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A RE-EXAMINATION OF RAHN’S DATA ON THE NUMBER 
OF GENES IN BACTERIA 


D. J. Finney 
University of Oxford 
L. MARTIN 


Institut Agronomique de Gembloux et Université de Bruxelles 


1. INTRODUCTION 


Fone to current cytogenetic theory, each chromosome in every 
cell of a living organism comprises a number of genes. Inter alia, 
these genes have the following properties:— 

(a) They are the carriers and transmitters of hereditary characters; 

(b) They control the specific enzymatic activity of the cell; 

(ec) They play an important part in the initiation of cell division. 
This paper is concerned only with the third of these properties, and 
with the light it can throw on the numbers of genes in a cell. The 
natural process of cell division must be preceded by duplication of the 
constituent chromosomes, which in turn requires the duplication of 
every gene and is the culmination of an elaborate synthetic process. 
Luria (1950) has pointed out that ‘There is something peculiar to 
homologous replication that sets it aside from other types of synthetic 
reactions. The replication of specific biological units must involve the 
building of complex specific molecules or molecular aggregates, the only 
permissible limitation to identity of model and replica being the pro- 
duction of “mutated” structures—the production, that is, of modified 
elements replicated in the modified form’. Claude (1949) has empha- 
sized the difficulties in this ‘template theory, elaborated to account for 
the supposedly unique process of gene reproduction’, which ‘assumes 
that the gene or the chromosome serves as a mold for the systematic 
apposition of a new substance’. He goes on to state that ‘in the light 
of biochemical processes already known, it is conceivable that the 
duplication of essential and characteristic cell substances is the end 
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result of a series of rigidly ordered chains of reactions, the final product 
in turn taking part at some point, and thereby directing the specificity, 
of the same or other, interlocked, biochemical cycles.’ 

In contrast to Rashevsky’s (1948) explanation of mitosis with the 
aid of a physico-chemical model of cytoplasmic behaviour, Rahn (1931) 
formulated a statistical model of the division of microorganisms in terms 
of gene duplication. Though his model was over-simplified, especially 
in that it related only to the nucleus and ignored cytoplasmic activity, 
it is important as indicating the type of phenomena to be expected if 
gene duplication is the chief controlling factor, and it ought then to be 
of value at least as an approximation to a more complex truth. Rahn’s 
discussion used the following assumptions about the nature of the genes 
in unicellular organisms:— 

(i) Each cell contains the same number, g, of genes; 

(ii) Each gene has a chemical structure that distinguishes it from 
other genes in the same cell and that it shares with genes at the corre- 
sponding chromosomal loci in other cells (genes corresponding in this 
manner will be referred to as genes of the same type); 

(iii) Gene duplication is a chemical process, and is subject to 
chemical laws; by analogy with the law of mass action, the rate of gene 
duplication at any instant is supposed proportional to the number of 
genes of the same type that have not yet doubled; 

(iv) The duplication of any gene takes place independently of that 
of all others in the same cell; 

(v) Cell mitosis takes place as soon as all genes in the cell have 
doubled; 

(vi) The probability that the gene at a particular locus doubles in 

any specified unit time interval is the same for all gene types. 
Of these assumptions, (i) would be generally accepted as a characteristic 
of the normal cell nucleus, (ii)-(v) are the basis of Rahn’s model, and 
(vi) is a simplifying assumption introduced (in the absence of any 
indications to the contrary) as an aid to the obtaining of mathematically 
tractable formulae. 


2. RAHN’S DISTRIBUTION OF TIMES 


For a specified cell, the probability that duplication at any locus has 
taken place within a time ¢ may be defined to be p(t). Assumption 
(iii) above shows that 


pt) =1-q, (1) 


where q is the probability that duplication does not take place in unit 
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time; by (vi), q is a constant for all loci. By (iv) and (v), the prob- 
ability that mitosis occurs within a time ¢, P(é), is the probability that 
the g genes of a cell have all doubled in that time: 


= (1 (2) 


For any possible values of g and q, P(t) increases steadily from zero at ee 
t = 0 to unity as a limit. Since P(é) is the cumulative frequency func- 
tion for mitosis, the probability of mitosis in time dé is f(t) dt, where 


i) = PO}. 3) 


It is convenient to replace the parameter g by a new parameter a, 
where 


q=e". (4) 
Hence 


= age "(1 (5) 


The probability of mitosis in any time interval is thus completely 
specified by two parameters, g and a. The genetic parameter, g, must 
be constant for one species of organism, whatever the conditions of 
experiment; the kinetic parameter, a, may vary from one set of observa- 
tions to another. 


3. METHODS OF ESTIMATING THE PARAMETERS 
Equation (2) may be written 
—at = log (1 — P’”’). (6) 


The ratio of the times required for any two specified values of P to be 
attained is independent of a, and could be the basis of an estimation of 
g in the manner proposed Fy Rahn. Using the agar hanging block 
technique, he had made counts of the times required for fission of 
individual bacteria in a culture of B. aerogenes. From this empirical 
assessment of the distribution P(é) he read the times that elapsed before 
10% and before 90% of the bacteria had divided, and estimated g from 
the ratio of these times with the aid of equation (6). Unfortunately, 
the ratio is not very sensitive to variations in g; Rahn was forced to 
conclude that ‘hence it seems rather hopeless to determine the number 
of gene-types by this method if the number is high, considering also the 
inaccuracy in determining experimentally the order of growth’. 

Ingenious though this method of estimation is, more sophisticated 
statistical techniques might be expected to yield better estimates of the 
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two parameters and also to allow tests of the adequacy of equation (2). 
In particular, the method of maximum likelihood is known to give 
estimates of maximal precision for large samples, under conditions that 
need not be specified heré (Fisher, 1922; Cramér, 1946). Even for small 
samples, the estimates are likely to be more precise than any based 
upon equation (6), because they make full use of the data instead of 
employing only certain percentiles. If the recorded times to fission for 
n cells are t; (¢ = 1, 2, --- , n), the likelihood function, L, is the product 
of expressions for f(t) and therefore 


log L = nloga+nlogg—a 


(2) 


The maximum likelihood estimates, g and a, are the solutions of the 
equations 


(8) 
of 
An iterative process for solving equations (8) may be set up (as for 


example by Finney, 1949; 1951, Appendix II), but this seems to con- 
verge rather slowly. An alternative is to write the equations as 


n 
9 > log u; 
and 
(9) 
=1+ 
where 
u,=1—e*". (10) 


The right-hand sides of equations (9) may be evaluated for two trial 
values of a, and a new trial value obtained by interpolation so as to 
make the difference between these zero. Further interpolations are 
made until both right-hand sides agree in giving the same value, which 
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is g. If the second differential coefficients of log L with respect to g 
and « are evaluated for the particular values g and a, and are written 
as a matrix 

(a logL a” log 
ag ag da 


Vlog L log L 
| ag da da” 


the elements of the inverse of this matrix, with signs changed, are, for 
large n, asymptotically equal to the variances and covariance of g and 
a (Fisher, 1922). Whether these quantities are reasonable approxima- 
tions to the variances and covariance when 7 is small is unknown, but 
the corresponding assumption is often made in other applications of 
maximum likelihood estimation and no alternative method of obtaining 
any idea of the variances is at present known. On this supposition 


VQ) = E +(g - 1) {fa} ] 


V@®=5+D, (11) 
and 
where 


4, ANALYSIS OF DATA 


Rahn (1931) reported the results of nine sets of observations on 
B. aerogenes, these being recorded as numbers of fissions occurring in 
successive five-minute time intervals. Professor Rahn has very kindly 
made available to us two additional sets of observations, those of 
February 16 and 23, which formed part of the same series but were 
omitted from his paper because they were small and showed slower rates 
of fission; there appears to be no reason fur excluding them from the 
comprehensive analysis described here. The complete data are shown 
in Table 1. Estimates of g and a, obtained from each series by applica- 
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TABLE i 


NUMBERS OF FISSIONS OF B. aerogenes 


OBSERVED IN SUCCESSIVE FIVE MINUTE TIME INTERVALS 


Time interval 
Experiment | 2 3 5 6 7 s 9 10 | 11 | 12 | 13 | 14 |>14| Total 

Feb. 16 2 0 4: 4 5 3 5 3 1 27 
Feb. 17 2 13 12 7 7 2 1 44 
Feb. 23 2 6 10 12 1l 5 1 0 1 48 
Feb. 24 2 3 4 20 16 ll 2 1 0 1 60 
March 2 1 0 3 19 22 21 10 5 2 1 84 
March 3 1 3 12 20 38 19 18 8 4 1 1 0 1} 126 
March 6 4 12 30 19 15 1 0 2 1 84 
March 10 6 10 16 26 18 9 3 1 0 0 0 4 93 
March 17 6 23 43 30 9 0 1 112 
Nov. 12 13 | 60| 24 3 100 
Nov. 14 5| 13 4 5 2 1 30 

Total 1 7 61 | 172 | 158 | 164 | 107 68 38 | 19 5 2 1 5 | 808 


tion of equations (9), are summarized in Table 2, the standard errors 
quoted being the square roots of the variances given by equations (11). 
In the calculations, fissions recorded as having occurred in time in- 


tervals nos. 2, 3, 4, --- 


TABLE 2 


were assumed to have occurred after times 73, 


SUMMARY OF MAXIMUM LIKELIHOOD ESTIMATES OF PARAMETERS 


Rate of 
No. of gene [Correlation 
Experiment | cells, | No. of genes, g Kinetic duplication} coefficient 
n parameter, @ per 5 between 
minutes | and 
Feb. 16 27 22.44 10.3/ 0.473 +0.071| 0.377 0.909 
Feb. 17 44 534 383 0.877 + 0.105} 0.584 0.978 
Feb. 23 48 244 + 139 0.708 + 0.077| 0.507 0.967 
Feb. 24 60 63.0 + 23.4] 0.652 +0.061| 0.479 0.937 
March 2 84 29.44 0.652 +0.052| 0.479 0.910 
March 3 126 37.5+ 9.0] 0.600 + 0.041 0.451 0.929 
March 6 84 89.9 + 32.6] 0.833 + 0.070| 0.565 0.953 
March 10 93 31.14 8.8] 0.593 +0.050| 0.447 0.931 
March 17 112 59.1+ 16.5] 0.991 +0.071| 0.629 0.941 
Nov. 12 100 | 1,010 + 493 1.594 +0.118| 0.797 0.979 
Nov. 14 30 91.5 + 60.1/ 0.997 +0.151| 0.632 0.961 
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123, 173, --- minutes; the grouping appeared to be sufficiently fine to 
make more elaborate treatment unnecessary. The five times beyond 
the fourteenth interval were assumed to be 723 minutes. 

For each experiment, expected frequencies in each time interval may 
be calculated by inserting g and @ in equation (2). Table 3 summarizes 
these calculations for the experiment of March 3. The data conform 


TABLE 3 


OBSERVED AND EXPECTED FREQUENCIES FOR THE EXPERIMENT OF MARCH 3 
(g = 37.45, @ = 0.6001) 


Time interval Proportion doubled Number of fissions in interval 
qj) at end of interval 
[P(57)] Observed Expected 
3 0.0012 1 0.2 
5 0.1478 12 15.0 
6 0.3546 20 26.1 
7 0.5681 38 26.9 
8 0.7340 19 20.9 
9 0.8442 18 13.9 
10 0.9113 8 8.4 
a4 0.9503 4 4.9 
12 0.9724 1 P| 
13 0.9848 1 1.6 
14 0.9916 of? 
>14 1.0000 1 1.0 
Totals 126 126.0 


well to expectation, the value of x” being 9.9 with 6 degrees of freedom 
(nine classes were used, expected and observed totals were made to 
agree, and two parameters were estimated). The other series would 
show similar good agreement if they were tested in this way. 


5. DISCUSSION OF DATA 


Chief interest attaches to comparisons between the eleven experi- 
ments. As emphasized above, the true value of the genetic parameter 
must be the same in all experiments. At first sight, the finding of 
estimates of g ranging from 22 to 1,010 seems to contradict this so 
strongly that, in spite of the large standard errors, the hypothesis and 
assumptions on which equation (2) rests must be wrong. Surely any 
test of significance for differences in the g column of Table 2 would show 
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them to be far too great for explanation in terms of sampling variation! 
The answer, however, is not quite as simple as at first appears. 

In the distribution of times represented by equation (2) or equation 
(5), the roles of g and a are biologically distinct but mathematically 
similar: each has an important influence on the location of the distribu- 
tion. Evaluation of the expected mean time gives 


(13) 
1 1 1 1 
+3) 
which, for large g, may be written 
E(t) ~ (y + log. g)/a, (14) 
where y is Euler’s constant, 0.5772 - - - ; even for g = 10, the error 


in replacing (13) by (14) is less than 2%, and it rapidly decreases for 
larger g, so that for most practical’ purposes (14) may be used. Conse- 
quently, pairs of estimates of g, a which keep the right-hand side of 
(14) constant will be equally satisfactory as judged by the first moment 
of the distribution; only the dispersion in the times can discriminate 
between effects of changes in g and of changes in a. This property of 
the distribution is evident during the process of maximum likelihood 
estimation: a considerable change in the value of one parameter can 
be made without seriously departing from the conditions expressed by 
equations (9), provided only that it is compensated by a change of the 
other parameter in the same direction. In fact, the covariance of the 
estimates is large relative to the variances; a correlation coefficient may 
be defined as 


_ CG, a) 
{V(g) 
and is found to exceed 0.9 for every experiment (Table 2). Moreover, 
the eleven pairs of estimates in Table 3 show the same tendency—a 
strong positive correlation. 
At our invitation, Mr. R. M. Jones has made a further study of the 
properties of the distribution represented by equation (5). He has 
shown that the cumulants have a simple asymptotic form when g is 
large: 


(15) 


6,& 


s=1 
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k, being given by (14). Even when g is only of the order of 30, these 
expressions are very near to the truth, and thus all cumulants except 
the first are nearly independent of g. In particular, the variance of the 
distribution of times is 

V(t) ~ 2/60’. (17) 


The first two moments of a sample of n observations might therefore 
be used for giving ap , go , estimates of the two parameters according to 
the ‘method of moments’. These would be defined by 


n-1 


6 
(18) 


Calculation of a and g, is easier than that of a and g. In nine of the 
experiments, a) is greater than @, and, as is to be expected, g, then 
exceeds g. The values of a range from 0.605 to 1.879, and those of go 
from 28.2 to 3,630. The method of moments appears to be appreciably 
less efficient than the method of maximum likelihood, however, since 
the sampling variance of a)» may be shown to be approximately 


2 
= (19) 


and the square root of this always exceeds the standard error ot @ in 
Table 2. 

In view of the correlation between. g and a, it is natural to inquire 
whether, despite their seeming discrepancies, the experiments may not 
agree satisfactorily with predictions based upon a common value of g. 
The value of a need not be the same for all, and the ideal estimation 
procedure would estimate g and eleven separate values of a from all the 
data. This would be exceedingly laborious, and, as an alternative, the 
more restricted theory that a also must be constant has been tried. 
Estimates of g and a found by applying equations (9) to the total line 
of Table 1 are 


g = 25.95 + 2.39 


& = 0.6139 + 0.0177, 


eee 

Qo (t) 

1 
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with 
Vs = 0 .924 


These estimates are only a little greater than the smallest recorded in 


Table 2; the explanation lies in the fact that the small values in that 
Table have smaller variances than the large, and therefore exert a 
dominant influence when the data are pooled. As confirmation that the 
smallness of g was not a consequence of serious bias in the maximum 
likelihood process, estimates by the method of moments were also 
formed, and are 


Jo = 34.2, 


a = 0.657 +*0.024. 


These are appreciably larger than g and a, but of the same order of 
magnitude; judging by a , they are less precise estimates. Table 4 is 
analogous to Table 3 and shows the observed and expected frequencies 


TABLE 4 


OBSERVED AND EXPECTED FREQUENCIES FOR THE POOLED DATA 
(9 = 25.95, & = 0.6139) 


Time interval Proportion doubled Number of fissions in interval 
(j) at end of interval 
[P(5j)} Observed Expected 
2 0.0001 
3 0.0113 7 9.0 
4 0.0975 61 69.6 
5 0.2911 172 156.4 
6 0.5165 158 182.1 
7 0.7008 164 148.9 
8 0.8255 107 100.8 
9 0.9016 68 61.5 
10 0.9455 38 35.5 
11 0.9701 19 19.9 
12 0.9837 5 11.0 
13 0.9912 2 6.1 
14 0.9952 1 
>14 1.0000 . 5 3.9 
Totals _ 808 808.0 


for the pooled data; the value of x’ is 14.9 with 9 degrees of freedom. 
This scarcely exceeds the 10% significance level (14.7), and does not 
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suggest any serious disagreement between theory and observation. If 
the unnecessary restriction on a were dropped, and each experiment 
were allowed its own kinetic parameter, the agreement would presum- 
ably be improved; the improvement might not be very great, because 
of the correlation between estimates. 


6. SUMMARY AND CONCLUSIONS 


Rahn suggested a theory of cell multiplication, in terms of the 
duplication of all genes in a cell, according to which the frequency of 
mitosis depends upon the number of genes per cell and a parameter 
representing speed of gene duplication (for simplicity, supposed constant 
for all genes). The distribution of times to fission that results from this 
theory has been discussed; a maximum likelihood procedure for the 
estimation of the parameters from observational data has been estab- 
lished. 

On the assumption that Rahn’s data for B. aerogenes are random 
samples of the times that elapse before fission of individual bacteria, 
estimates of the parameters have been formed from his eleven experi- 
ments. All agree satisfactorily with the theory; though the separate 
experiments show large differences in their estimates, when all the data 
are combined they agree well with the theory if the number of genes 
per cell is about 26 and the rate of gene duplication is about 46% per 
five minutes. The theory actually tested has been unduly restrictive 
in demanding that the kinetic parameter be the same for all gene-types 
and on all occasions. Lesser requirements would obviously still be in 
accord with Rahn’s physicochemical hypothesis, but the data are in- 
sufficient for investigation of whether q is really constant on all occasions 
and for all genes. Though statistical tests disclose no significant de- 
parture from the predictions based on equation (2), the estimate of the 
number of genes is so remarkably small as itself to cause suspicion. In 
correspondence with one of us, Professor Rahn commented ‘I am as 
surprised about the small number of genes that you found as I was 
about my own results. To the uninitiated, a bacterium may seem a 
very simple organism, but to the bacteriologist it seems almost un- 
believable that 30 genes could control all the many properties of such 
a cell. There is not only the ability to make cell wall and protoplasm, 
but the ability to decompose carbohydrates to various acids, carbon 
dioxide, and hydrogen, and above all to construct all of the catalysts 
needed for growth, e.g. the mechanism to make various amino acids 
from glucose and ammonia and to make proteins and nucleoproteins out 
of these amino acids. It must be that each gene controls a number of 
synthetic properties.’ The only alternative to the view expressed by 
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Professor Rahn’s last sentence appears to be that the theory expressed 
by equation (2) is seriously wrong in the importance it attaches to the 
number of genes as a factor determining rate of fission, but that the 
data analysed happen to agree with it as well as with the unknown 
correct theory. 

The problem of estimating the number of genes in a micro-organism 
is essentially biometric in character (cf. Cochran, 1950; Martin, 1949), 
involving the formulation of a mathematical model of certain biological 
phenomena and the submission of this model to experimental test. 
Unfortunately, the mathematical model appropriate to Rahn’s theory 
is such that a reasonably sensitive test of its adequacy is impossible 
without either vastly more data than he obtained or experimentation 
of an entirely different character. This conclusion, the clearest that 
emerges from the present paper, can give little satisfaction to statisti- 
cians or to experimenters, and, in the face of it, the numerical estimates 
must be accepted with reserve. 
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RECTANGULAR LATTICES AND PARTIALLY BALANCED 
INCOMPLETE BLOCK DESIGNS 


K. R. Narr 
Forest Research Institute, Dehra Dun, India 


 aecmamonapnag (1946, 1947, 1949) developed a set of incomplete block 
designs in which the number of varieties (or treatments) is the 
product of two consecutive integers and the number of replications of 
each variety is either 2 or 3 and their multiples. He called them simple 
and triple rectangular lattices respectively. 

It is the purpose of this paper to show that the simple rectangular 
lattice belongs to the category of partially balanced incomplete block 
designs developed by Bose and Nair (1939) and later modified by Nair 
and Rao (1942). It will also be shown that the triple rectangular lattice 
does not, in general, satisfy the combinatorial conditions required of a 
p.b.i.b. design. 

A brief reference to these conditions had been made in this journal 
by Bose (1947). They may be stated as follows:— 

(1) There are v varieties (or treatments) each replicated in r plots. 

The vr plots are arranged on the ground in b blocks each containing 

k plots (k < v). No treatment will occur more than once in a block. 

(2) With respect to any treatment 6 the remaining ones can be 

divided into m groups of n; , n2, -*- , 2», each, so that each treatment 

of the i-th group appears with @ in \,; blocks. The treatments of 
. the i-th group are said to be i-th associates of 0. - 

(3) The relation of association is symmetrical. Thus if the treat- 

ment @ is an i-th associate of ¢, then ¢ is an i-th associate of 6. 

Further, if pj, is the number of treatments which are both j-th 

associates of 6 and k*-th associates of ¢, then the value of p}, is inde- 

pendent of the pair of i-th associates @ and ¢ with which we start. 

The following relations exist:— 


*This k is not to be confused with the number of plots per block which is also denoted by k. 
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vr = bk (1) 
3 nar; = r(k — 1) (2) 
Dix Dii (3) 


Dis =n; —1 or n; 
(4) 


according asi = j or tj 


= = Mi; (5) 


Originally, Bose and Nair (1939) had imposed the condition that 
A.» *** » Am Were unequal, but Nair and Rao (1942) found that this 
condition was not necessary. 

The treatment differences are estimated with m different accuracies. 
The ‘balanced incomplete block design’ can be treated as a special 
case of p.b.i.b. designs, with m = 1. The square and cubic lattices for 
p and p° varieties respectively, in blocks of p plots, are special cases 
of p.b.i.b. designs, with m = 2 and 3. Many illustrative examples of 
these designs for which m is 2 or 3 have been given by Bose and Nair 
(1939). It was thought at the time that for practically useful designs 
m should not exceed 3. ; 

The simple rectangular lattice for p(p — 1) varieties is a fine example 
of a p.b.i.b. design for which m = 4 when p > 4, and m = 3 when p = 3. 
Since only two replications are required for this design its practical 
usefulness is obvious. 

Following Harshbarger’s method, each of the p(p — 1) varieties 
will be identified by a couple of digits x and y (x ¥ y) each taking the 
values 1, 2, --- p. All varieties having the same value for x will form 
the z-th block of the first replication and all varieties having the same 
value for y will form the y-th block of the second replication. 

For any particular variety (cy) there are 2(p — 2) varieties ap- 
pearing with it in one or the other of the two blocks in which it is repli- 
cated. Each of these 2(p — 2) varieties has one digit in common with 
variety (xy). They will be called Ist associates of (ry) with A, = 1. 
Hence n, = 2(p — 2). 

The remaining p’ — 3p + 8 varieties can be grouped into three (if 
p > 4) according as they have 0, 1 or 2 digits in common with (ay). 
They may be called 2nd, 3rd, and 4th associates of (xy) with A. =. 0, 
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As = 0, Ay = Oand n, = (p — 2)(p — 3), n, = 2(p — 2) and nm = 1 
respectively. 


Values of pj, can be easily worked out. Their values and those of 
the parameters obtained above are given below for p= 4. 


v=p(p—1), k=(p—1), r= 2, b = 2p 
A = 1 = 0 | =0 


m = 2p — 2), m = (p — 2)(p— 3), ny = Ap— 2), m = 1. 


[ (p — 3) (p — 3) 1 0 
(p—3) (p-—3) 0 
@-3)  @-3) 1 
0 0 1 0) 
[2 2(p — 4) 2 0) 
2p-4) (-4(p-5) Ap-4) O 
aah ae" 2(p — 4) 2 0 
[0 1 0 0) 
. 4 (p — 3) (p — 3) 1) 
(P—3) (p—3) 1 
0 0 
( 0 0 0) 
; 0 (p — 2)(p — 3) 0 0 
~~ 2(p — 2) 0 0 0 
| 0 0 0 0 


i 
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When p = 3, n, becomes zero so that m reduces to 3. The 3rd and 


4th associates for the case p > 4 may be called 2nd and 3rd associates 
respectively for the case p = 3. The parameters are 


ink veh 
A, = 1 A, = 0 A, = 0 


n, = 2, N, = 2, nm = 1 


010 0 2 0 
Diz = 1 0 1 ’ Pa = ef 0 ’ Pix = 2 0 0 


The statistical analysis appropriate for p.b.i.b. designs was discussed 
by Bose and Nair (1939) taking into account only intra-block informa- 
tion. Methods for using inter-block in addition to intra-block informa- 
tion have been given by Nair (1944) and Rao (1947). 

For the simple rectangular lattice, Harshbarger directly derived 
expressions for the various sums of squares in the analysis of variance 
and for variances of treatment differences. Since this lattice has now 
been shown to be only a special case of p.b.i.b. designs these expressions 
can be obtained by substitution in the general expressions for the p.b.i.b. 
designs. 


Triple rectangular lattice. I had naturally expected that, like the 
simple rectangular lattice, the triple rectangular lattice discussed by 
Harshbarger (1949) will also be a special case of p.b.i.b. designs, with 
m taking the value 7 in general. But this expectation has not come 
true except in the cases v = 2 X 3 and v = 3 X 4 when they can be 
shown to be special cases of p.b.i.b. designs with m equal to 2 and 3 
respectively. 

In the triple rectangular lattice each of the p(p — 1) varieties is 
identified by a set of 3 digits ryz (x * y ¥ 2), each taking the values 
1 to p. All varieties having the same value for z will form the z-th 
block of the first replication, all varieties having the same value for y 
will form the y-th block of the second replication and all varieties having 
the same value for z will form the z-th block of the third replication. _ 

For any particular variety (xyz) there are 3(p — 2) varieties ap- 
pearing with it in one or the other of the three blocks in which it is 


4 
ae 
4 
oa 
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replicated. These 3(p — 2) varieties can be grouped into three (if 
p > 5) according as they have 1, 2 or 3 digits in common with (xyz). 
If p = 4, the number of groups will be two according as the number of 
common digits is 2 or 3. If p = 3, there is only one group, as all the 
3 digits are common for every variety. 

The remaining p> — 4p + 5 varieties can be grouped into four (if 
p = 6) according as they have 0, 1, 2 or 3 digits in common with (xyz). 
If p = 5, there can be no two varieties with no digits in common and 
hence the number of groups reduces to three. Similarly when p = 4, 
there will be only two groups, namely those with 2 or 3 digits in com- 
mon. Lastly, when p = 3 there is only one group as all varieties have 
all the 3 digits alike. 

In general, therefore, when p > 6, with respect to any variety (xyz), 
the remaining varieties apparently fall in seven groups. 

Let us take the smallest triple rectangular lattice, namely, that for 
which p = 3. The six varieties may be written in the form 


x 123 132 
213 XxX 231 
312 321 x 


It can easily be seen that this is a p.b.i.b. design with following 
parameters:— 


v = 6, k = 2, r= 3, b = 9. 


=2 


0 2 3 0 
Dit Dir 
2 0 


Let us next take the case p = 4. The 12 varieties may be written 


\ 
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in the form 
x 123 134 142 
214 x 231 243 
312 324 x 341 
413 421 432 x 


We notice that with respect to any variety all the varieties which 
appear with it in the same block have 2 digits alike; and among the 
varieties which do not so appear some have 2 digits alike and others 
have 3 digits alike. If we call these groups the 1st, 2nd and 3rd asso- 
ciates respectively, the above triple rectangular lattice is a p.b.i.b. 
design with following parameters:— 


vy = 12, k = 3, r= 3, b = 12. 
= 1 = 0 4s = 0 
n, = 6 %=3 = 2 
221 40 2 3 3 0 
0 1); Pr= 10 2 Pr = |3 0 0 
1 10 20 0 001 


Harshbarger (1949) gives a numerical illustration of this lattice. On 
p. 11 of his paper five different standard errors are worked out of which 
three are said to be for differences between the means of two varieties 
occurring together in the same block and 


(1) having one digit alike in the subscripts; 
(2) having two digits alike in the subscripts; 
(3) having three digits alike in the subscripts. 


@ 
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It is clear from the discussion above that groups (1) and (3) are 
non-existent. 

The remaining two standard errors worked out by Harshbarger are 
said to be for differences between the means of two varieties not occurring 
together in the same block and 


(1) having one digit alike in the subscripts; 
(2) having two digits alike in the subscripts. 


Of these, (1) is non-existent. Furthermore among varieties not oc- 
curring together in the same block there are pairs having three digits 
alike which have not been taken into account by Harshbarger in the 
evaluation of standard errors. 

We will now examine the triple rectangular lattice for p = 5 and 
show that it does not satisfy the condition for a p.b.i.b. design. The 
20 varieties may be written in the form 


Xx 124 132 145 153 
214 x 235 243 251 
312 325 xX 341 354 
415 423 431 x 452 
513 521 534 542 x 


Take variety 124. The remaining 19 varieties can be grouped into 
five classes as follows:— 
(a) varieties appearing with 124 in the same block and having either 
1, 2, or 3 digits alike; 
(b) varieties not appearing with 124 in the same block and having 
either 1 or 2 digits alike. 
The varieties falling in the five groups are:— 


| 
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Group: a(1) a(2) a(3) b(1) b(2) 
153 132 214 235 243 

325 145 513 251 

354 423 312 

534 521 341 

415 

431 

452 

542 


If this lattice was a p.b.i.b. design these five groups would have re- 
spectively formed the 1st, 2nd, 3rd, 4th and 5th associate classes of 
variety 124 and we would have had A, = A, = Az = 1,4 = As = O 


and n, = n, = 4,n, = 1,m = 2, n; = 8. 


Assuming that it is a p.b.i.b. design let us examine the values of the 
parameters p;,. Take, for instance, variety 534 which is a 1st associate 


of 124. The varieties belonging to each associate class of 534 are: 


Associate 
No. 1 2 3 4 5 
521 513 354 251 145 
132 542 312 153 
124 235 243 
214 431 325 
341 
415 
423 
452 
We can now evaluate the values of p;, . They are 

(0 0 1 0 2) 

200 0 2 

000 0 

0200 0 

(0 2 0 2 4) 
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We find that p}, ¥ pi; for some values of j and k. Hence condition 
(3) of a p.b.i.b. design is not fulfilled. By successively taking a variety 
each from the 2nd, 3rd, 4th and 5th associate class of 124 and trying 
to evaluate pj. , Pix , Dix and p;, it can be shown that condition (5) of 
a p.b.i.b. design also breaks down for the 4 X 5 triple rectangular lattice. 

When we come to 5 X 6 lattice further asymmetries appear in the 
design. For instance, consider the 30 varieties numbered as shown 
below: 


x 123 134 145 156 162 
216 x 231 243 254 265 
315 326 x 341 352 364 
412 425 436 x 451 463 
513 524 532 546 x 561 
614 621 635 642 653 x 


It can be seen, for instance, that the number of varieties having 3 
digits alike with variety (xyz) is either 0, 1 or 2 and not constant. 
Examples of the 3 types are (xyz) = 243, 123 and 254 respectively. 

Take another 5 X 6 lattice where the z-digits are arranged in a 
different way as is shown in the next diagram. Here the number of 
varieties having 3 digits alike with variety (xyz) is either 0 orl . Examples 
are (xyz) = 241 and 126 respectively. 

Such asymmetries also appear in lattices for p = 8 and 10 listed 
by Cochran and Cox (1950). 

We had seen that the number of varieties having 3 digits alike with 
variety (xyz) is 1 when p = 5. This is true also for larger odd values 
of p. Hence the number of groups of varieties with respect to any 
variety (xyz) is never 7 when p is an odd integer. 
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x 126 132 143 154 165 
213 x 235 241 256 264 
314 325 x 346 351 362 
415 423 436 x 452 461 
516 521 534 542 x 563 
612 624 631 645 653 x 


The fact that the triple rectangular lattices for p > 5 are not 
p-b.i.b. designs seems to show that there is a more general class of 
partially balanced designs with less stringent combinatorial conditions 
than those imposed in the p.b.i.b. designs developed by Bose, Nair 
and Rao. The combinatorial conditions required for this more general 
class of designs have not been looked into. 
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NOTES ON THE ESTIMATION OF GROSS AND NET 
REPRODUCTION RATES BY METHODS OF 
STATISTICAL SAMPLING 


J. C. Koop 
Directorate of Labour, Rangoon, Burma 


1, INTRODUCTION 


” gay AND NET reproduction rates are usually computed from data 
on births and deaths relating to the total population. But if the 
data are faulty the degree to which such figures are correct cannot be 
assessed. The purpose of these notes is to point out that these rates 
may also be estimated from data, possible to obtain, by methods of 
statistical sampling which also supply methods of estimating their 
accuracy. In the case of the latter index, supplementary actuarial in- 
formation in the form of female survival rates, l, , for specific ages are 
necessary. 

If, at a given time period, in any geographical region the entire 
population of women is considered, the gross reproduction rate G, ac- 
cording to Kuczynski, is simply the sum of the specific fertility rates 
for female births p, , where the subscript 2 denotes, the age-group 
(« — 4, x + 3) of the mother. The net reproduction rate R invented 
by Béckh (D.V.G., 1947), will be the sum of the products of the age- 
specific fertility rates and each corresponding female survival rate. 
We need only be concerned with those groups included in the repro- 
ductive period which is usually between 15 and 50 years. In symbols 


49.5 49.5 
= > p. and R= >> lp. 
z=15.5 z=15.5 


These formulae are frequently used by demographers. In this problem, 
each p, may be regarded as a population parameter to be estimated 
by sampling, and each /, as a known constant. 


2. SAMPLING FROM A FEMALE POPULATION STRATIFIED IN AGE-GROUPS 


Suppose in a population, all females between 15 and 50 years are 
located, listed and stratified in yearly age-groups. Such an event can 
be realized when a census is being taken. Let N, be the number of 
females in the age-group (x — 3,x + 3). Treat each female as a samp- 
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ling unit and assign to each an equal probability of being selected. If 
random samples of n, females are selected from each age-group, and 
m, are found to be mothers of female infants under one year in the 
age-group (x — 3, x + 4) then an unbiassed estimate p! of p, will be 
given by p, = m,/n,. Hence unbiassed estimates G’ of G and R’ of R 
will be 


and R= 


Needless to say these estimates will be subject to sampling error. 
We now proceed to find, their variances which measure these errors. 
With each female, we may associate a characteristic random variable 
which assumes the value one, if the female in question is a mother of 
baby girl of age less than one year, and zero if otherwise. With this 
notion, and also remembering that variables not in the same age-group _ 
are mutually independent, but those within dependent, it can be easily 
shown that the variances of G’ and R’ are 


(1) 
and 
V(R’) = Non (2) 
Their estimates are: 
1 of 
respectively, where f, = n,/N, , the sampling fraction in each age- 


group. These estimates are unbiassed since their expected values are 
V(Q@’) and V(R’) respectively. If the number N, in each age-group is 
effectively large, f. may be neglected. 

To carry ‘this discussion further, if circumstances are such that a 
sample as large as m, = > 2*15.5 m- can only be selected, it can be 
shown that the variance of G’ will be a minimum if the age-structure of 


the sample is such that 


N, 
pAl — p.) N. V/pAl — pz) 


This of course is the analogue of the formula for quantitative variates 


(for large N,’s). (3) 


z =z 
| 
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given by Neyman (1934). To prove the above result, it is necessary 
to minimize V(G’) with respect to the n,’s subject to the restriction 
n= b Dene n,. This may be done by the method of Lagrange multi- 
pliers. If cost is a practical consideration, and if it costs ¢ units to in- 
vestigate every female sampled, and if monetary resources are limited 
to @ units, then a/c females can only be sampled. The age-structure of 
the sample for which V(G’) is a minimum, subject to the cost restriction, 
will be given by substituting a/c for n, in (3). Similarly under the 
same, size-, or cost-restrictions, n2 = >, n, the variance of R’ will be 
a minimum if the allocation of sampling units among age-groups is 
such that 


Ne 
N, — pz) 
(4) 
(for large N,’s). 


If the problem is considered from the point of view of finding the 
optimum age-structure of a sample, subject to the restriction that the 
variance shall be a fixed quantity, the solutions obtained in respect of 
either G or R are identical with those obtained above. 

We may now consider the implications of equations (1), (2), (3) 
and (4): 

(i) It must be admitted that there is no exact solution to the problem 
of determining sample size since the p,’s are themselves the quantities 
we are attempting to estimate at the time period in question. However, 
to obtain a workable solution, the p, values of a previous time period, 
not very far off, may be used. In the discussion which follows we shall 
bear this in mind. 

(ii) Consider the problem of sampling so as to have either the 
estimate of G or R with minimum variance. If an investigator wishes 
to estimate both these rates from the data of a single sample, a de- 
cision will have to be made before the sample units are allocated among 
the various age-groups, as to whether G’ or RP’ is required with minimum 
variance, for, an inspection of (3) and (4) shows that the age-structure 
of a sample of a given size which confers the former with a minimum 
variance does not’ do similarly to the latter and vice versa. The fol- 


1The only trivial case which is an exception to this statement is when /z = a constant for all age- 
groups under reference. This implies a state of affairs in which all females between 15 and 50 years 
have equal chances of survival. In the present condition of mankind this possibility could not be 
hoped for. 
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lowing question almost suggests itself: how may control be exercised 
over the estimate which is not given a minimum variance? We now 
proceed to answer this question. If the sample is designed to give R’ a 
minimum variance, then its age-structure will be as given in (4), so 
that in this situation, the variance of the estimate of @ will be <aing by 


where n is the size of the sample and the symbols W, , W. and W, 
(which have been introduced for the sake of brevity), may be identified 
in the order presented in the identity. On the other hand, had the age- 
structure of the sample been according as given by (3), the estimate 
of G would have had a minimum variance of 


1 N, 1 = 


N, 


whence the relative efficiency of this estimate would be given by 
= (W; — nW)/(W,W, — nW), (5a) 


and, in view of the fact that W — 0, for large N,’s, its limiting value 
will be given by 


(5b) 
1 
L — Py. 
Clearly E, cannot exceed E,, so that this limiting relative efficiency, 


which is independent of sample size, can only be attained in the case 
of samples from effectively large populations. We also have 


Vo — Von = (WiW2 — W2)/n 
N rez’ (6) 
(i each we ; is taken as 1). 


The expression W,W, — Ws is a function of p, 1 and N, and on the 
assumption that some prior knowledge of these quantities exists, suffi- 


N,-1 
| 


ESTIMATION OF REPRODUCTION RATES 159 


cient to obtain an approximation, it may be treated as a known constant. 
It follows that 6 = V_ — Ve, , the deviation from the minimum 
variance, is inversely proportional to the sample size n. Hence for any 
given 6, the sample size n can be determined approximately. Its 
accuracy of course will depend on how closely W,W, — W3 approxi- 
mates to its true value. We see that if these n sample units are allocated 
among the age-groups in the way indicated by (4), R’ may be expected 
to have a minimum variance, and simultaneously, G’ may also have a 
variance differing somewhat closely from its ideal or possible minimum 
variance. Hence some degree of control over @’ has been achieved. 
This approach has been adopted in view of the fact that demographers 
attach more importance to the net reproduction rate. With the alterna- 
tive approach, when G’ is given a minimum variance, the relative 
efficiency of the estimate of R, and the deviation of the variance of the 
estimate of R from its ideal minimum variance are given by 


Ep, = (Wi — nW,)/(W2W; — nW,), (7) 

and 
Ve — Vim = (W3W; — Wi)/n 

1 

(1 — p.)(L — De’) (8) 
) 
(i each V.-1% taken as 1 

where 


W,= > p:)/(N.z 1), 


W; 


N. 
2 


and n, W, and W; are as defined above. 

(iii) In the alternative problem, when the allowable amount of 
variance is fixed in advance, the investigator who wishes to estimate 
both G and R from the same data will again encounter the problem of 
controlling the estimate whose variance has been left unfixed. Suppose 
for example that the variance of R’ is fixed at v. Remembering that 
the age-structure of this sample will be given according to (4), we have 


w?-W,. (9) 
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Since the conditions which give rise to equation (6) are formally equiva- 
lent to those in the present problem, we have 


5 = (W,W, — W3)/n 


By eliminating » between (9) and the above equation, the relation 
between v and 6 is seen to be 


WwW; Wr 


And on the assumption that the expressions involving the W’s can be 
approximated from prior information, (10) may be used to make nu- 
merical adjustments between the deviation from the ideal minimum 
variance of G’, and, the “fixed” variance of R’. After such adjustments 
are made n can be determined by substituting the value of the “fixed” 
v in (9). Also since the formulae for the allocation of sampling units 
in both problems are identical, the relative efficiency of the estimate of 
G will be as given in equation (5a). Hence we see that some control 
over the unfixed variance can also be exercised. Similar formulae may 
be developed for the case when the variance of G’ is “‘fixed’’ and that 
of R’ unfixed. 

(iv) Lastly it may be said that the frame of this sampling system, 
in which the unit of sampling and the unit of analysis are identical, is 
perhaps not practicable for large populations of the order of millions. 
This is true. But for the case of small ethnic groups or minority popu- 
lations, living within compact geographical areas and in whose prospects 
of biological survival we may be interested, these methods may be 
suitable. These populations may number anything between five and 
fifty thousand. It is also for this reason that the above theory has been 
developed with a finite approach. 


(10) 


3. ILLUSTRATIVE EXAMPLE 


An example to illustrate the above theory on the optimum alloca- 
tion of sampling units in the estimation of the gross reproduction rate 
now follows. The data which is presented in columns (i), (ii) and (iii) 
of the table is drawn from a report on a recent census’ of the Eurasian 


2This census, which was sponsored by the Anglo-Burman Union, was carried out by the writer 
with the assistance of a group of honorary workers in 1949. The writer is grateful to H. A. Faria and 
W. E. Foster for their generous help in the work of enumeration. Due to the conditions of civil war 
around Rangoon, the census could not be completed. The Eurasian community numbered about 
6,200 at that time. The data on fertility rates is based on 81% of the total female population in the 
reproductive group. There was also one other defect. It was neglected to ask whether female daughters 
of age less than 12 months had died in the 12-month period prior to the census. For those age-groups 
where this defect is present, it is possible that fertility rates may tend to be underestimated, unless 
compensated by a syst tic tendency to select households containing mothers whose female infants 
survived the first year of life. However, I am of the opinion that this defect is not serious since infantile 
mortality among Eurasians is not high. 
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people in Rangoon, which so far has not been published. 
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(Since sur- 


vival rates of Eurasian women are not known, it will not be possible 
to give illustrative examples on the problems of sample size, when both 
the net and the gross reproduction rates are required to be estimated 
from the data of a single sample.) 


TABLE 
Age- | Number |Female infants} Age-specific V — pz) 
group of below fertility V — pz) ZV — ps) 
emales one year rate pz 
(i) (ii) (iii) (iv) (v) (vi) 
17-18 47 2 .042553 20185 .02925 
18-19 62 0 20185! .02925 
19-20 57 3 .052632 22330 .03236 
20-21 66 7 . 106061 30792 .04462 
21-22 60 2 .033333 17950 .02601 
22-23 57 5 .087719 28289 .04100 
23-24 52 4 .076923 26648 .03862 
24-25 46 4 .086957 28231 .04091 
25-26 51 5 .098039 29918 .04336 
26-27 48 8 . 166667 36268 .05259 
27-28 38 2 .052632 22330 .03236 
28-29 57 2 .035088 18400 .02666 
29-30 57 5 .087719 28289 .04100 
30-31 52 4 .076923 26669 .03865 
31-32 25 1 .040000 19596 .02840 
32-33 41 4 097591 29676 .04301 
33-34 30 2 .066667 24944 .03615 
34-35 35 a .028581 16667 .02415 
35-36 37 1 027027 . 16216 .02350 
36-37 36 2 055556 . 22906 .03319 
37-38 28 1 035714 . 18558 .02689 
38-33 31 0 18558 |. .02689 
39-40 34 0 18558) .02689 
40-41 Ad 2 .045455 20830 .03019 
41-42 21 0 20830 .03019 
42-43 26 0 20830 .03019 
43-44 20 0 20830 3 .03019 
44-45 23 0 20830 .03019 
45-46 39 0 . 20830 .03019 
46-47 18 1 .055556 . 22906 .03319 
G = 1.455393) 6.90049 1.00004 

1Token figure, obtained on the assumption that pis.s = piz.s = .042553 

2Token figures, obtained on the assumption that pso.5 = p3s.s = p37.5 = .035714 

3Token figures, obtained on the assumption that pu.s = pus = = = 0045455 
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In column (iv) of the table, the age-specific fertility rates of Eurasian 
women are given. If these rates are added up, a figure which will 
somewhat approximate to the true gross reproduction rate will be 
obtained. 

If column (iii) of the table is inspected, it will be seen that no 
births are shown against women in certain age-groups. It is difficult 
to decide whether fertility rates for each of these age-groups should be 
taken as zero, since it was neglected to ask whether female daughters 
of age less than 12 months had died in the 12-month period prior to 
the census. If the procedure of gathering data had been statistically 
random, then there would be complete justification in taking these rates 
as zero. However, as the data of an incomplete census cannot be taken 
to be equivalent to that of a random sample, there would be no such 
justification. If it is presumed that these rates are zero, and the opti- 
-mum age-structure of the required sample is determined theoretically 
by formula (3), then no sample units will be allocated to these groups. 
In the event of their being not really zero, the sampled gross repro- 
duction rate will be underestimated, and hence biassed. For this 
reason, these groups should also be sampled. And it would be ex- 
pedient to assign a token rate to each of these groups. In assigning a 
token rate to any group, in this problem, the rate for the nearest group 
junior to the group in question has been chosen. This procedure. is 
almost arbitrary, and is founded on the hope that each of the fertility 
rates chosen for these groups might be about equal in magnitude to its 
corresponding ‘‘missing” rate. 

If the sample units are drawn at random from each age-group ac- 
cording to the proportions prescribed in column (vi) of the table, and 
second, if the age-specific fertility rates used for the purpose of de- 
termining these proportions happen to be just the rates prevailing at 
the period of time when the sampling is carried out, then the variance 
of the estimate of G will be a minimum. In regard to this problem, 
even if a sample survey had been carried out immediately after the 
census, the second set of conditions would not have been realized, for 
one reason, because of the fact that ‘“‘missing” rates had been somewhat 
arbitrarily chosen. However, selecting the sample units in this way, is 
definitely better than selecting them without any method, and even if 
G’ does not attain minimum variance, the experimenter will have 
grounds for belief that its variance is somewhere near the minimum. 


4. SAMPLING FROM A UNIVERSE OF LARGE PRIMARY UNITS 


As stated above, if the population is of the order of millions and 
spread over a very wide area, clearly the above sampling procedure 


| 
4 
4 
j 


ESTIMATION OF REPRODUCTION RATES 163 


will not be practicable. Hence the following alternative procedure is 
suggested. As in most sample surveys, suppose the area over which 
the population is distributed is unambiguously divided up along natural 
or adminsitrative boundaries into K small-size areas. Since the total 
area is supposed to be very large, the number of these primary units 
will also be very large, and each will also contain a large number of 
persons. Following Mahalanobis, these primary units may be termed 
grids. It may be noted that in each grid we are only interested in 
women in the reproductive age-groups, and their female children under 
one year. For a moment we shall defer general consideration on the 
question of what may be their best relative sizes, sex or age-composition. 
The notation adopted here will be substantially the same as that above, 
except for the addition of the suffix 7 to identify grids. Let n,,; be the 
number of females in the age-group (x — 3, x + 4) of the ith grid, 
(¢ = 1, 2, --- K), who have between them m,,; female children under 
one year. In respect of any given age-group in the universe of grids, 
it would be prudent to take into account the possibility that n,; and 
m,, may be correlated. The age-specific fertility rate will be given by 


K K 
This is the population value. The problem here is firstly to estimate 
each of these p,’s so as to enable us to estimate G and R, and secondly 
to estimate the sampling errors of these estimates. 

Let each grid be assigned an equal probability of selection. Let a 
random sample of k grids be selected and observe in each, the variates 
n and m, which have been defined above. An estimate p!’ of p, will 
be given by 


k 


k 
p = Yim: / Doni = 


i=l i=1 


This estimate, which is a ratio of two means, though consistent, is 
biassed, but Hansen and Hurwitz (1943) are of the opinion that the 
bias “is usually negligible’. Hence the estimates of G and R are as 
follows: 


od 
zn, z Nz 


The variance of G’’ can be shown to be 
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Vim.) _ Cov(m. ,n.) , Vin.) 


The expression (K — k)/(K — 1) accounts for the fact that we are 
sampling a finite population. We have remarked above on the possi- 
bility that m,; and n,, may be correlated. The covariance term in 
(11) accounts for this possibility. The estimate of V(G@’’) will be seen 
to be 


(1 — k/K) 
Kk — 1) Lip: E (m.; — ™,) 


k 
mn, i=1 
Correspondingly similar expressions may also be written down for the 
variance of R”. They will only differ in appearance by the term [2 
under the summation sign. 

It may be noted that the expression for the variance of G”’ given 
in (11) is only an approximate one. Actually it under-estimates the 
true variance. But it will suffice if each of the variances under the 
summation sign are small relative to their corresponding squared means 
(Yates, 1949). In this connection, Finney (1949), alluding to infinite 
normal populations, when the two variates comprising the ratio are 
independent, states that 7; should be at least 40 times V(7,) so as to 
be able to construct, without danger of error, 95% confidence limits 
for p, °. Since we are not attempting to find confidence limits for G, 
this condition perhaps in a less stringent form, applies also to finite 
populations. Hence to estimate the standard error of G” using the 
expression (11), conditions governing its validity must first be satisfied. 
To do this the only alternative seems to be to demarcate the grids in 
a way such that the numbers found in each reproductive age-group, 
ranging from 15-16 to 49-50, be 

(i) as near as constant from grid to grid, and 

(ii) sufficiently large, subject only to the restriction that the in- 
vestigator shall not find the grid too unwieldly to handle, or from the 
point of view of the planning authority, too costly to operate. 

(i) and (ii) ensure that for any given age-group, V(n,,;) will be con- 
trolled at a low figure and 7, at a high figure. Hence in large scale 


31 have given the substance of his statement using my own notation. 
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sampling, when k& will obviously be large, the ratio V(n,)/nz will be 
reduced to a very small value. Also since m,; is possibly correlated 
with n,; , V(m,)/m: may be expected to be indirectly controlled at a 
small value. 

The above statements presuppose that prior knowledge, both on the 
geographical distribution and age composition of the population, is 
available. In all developed countries, census records will furnish in- 
formation bearing on these questions. However, because of population 
changes, by way of births and deaths, and shifts from one locality to 
another, which are almost certain to occur, such records may not be 
completely satisfactory for the purpose in hand. But if a survey is 
timed to occur soon after a national census, then it may be possible to 
carry out the above suggestions fairly satisfactorily. If the age compo- 
sition of the female population is not known, but only the total numbers 
inhabiting each locality are known, the grids may be demarcated such 


_that the sex ratio in each grid and the numbers included in each grid 


remain almost constant, and also such that condition (ii) is satisfied. 
More often than not, this may be the only practical course open to the 
investigator. By this procedure it is quite possible that the females 
belonging to various age-groups, may be more or less equally distributed 
over all the grids. 

In this system of sampling, the cost of a survey will depend partly on 
the number of grids selected, and partly on the physical size and popu- 
lation size of each of the grids included in the sample. Hence if monetary 
resources are fixed at a certain limit, and if we have some empirical 
knowledge of average costs per grid, the number of grids to be selected 
will be determined by this limit. 

In closing this section it may be stated that alternative systems of 
sampling are certainly not lacking. Also, what has been suggested 
above may not be the best in all circumstances. The writer does not 
propose discussing these problems since the object here is merely to 
point out that gross and net reproduction rates may be estimated by 
methods of statistical sampling. Different systems of sampling have 
been discussed by various writers in the United States, Britain and 
India, and any of these systems may be used. 


5. NOTE ON pz 


A point arising out of this discussion concerns the interval estima- 
tion of p, , the age-specific fertility rate for female births. Suppose that 
large samples are being selected from very large populations. If p, is 
being estimated by the system of sampling outlined in section 2, then 
its confidence limits, corresponding to the confidence coefficients 0.95 


: 


166 BIOMETRICS, JUNE 1951 


or 0.99, for certain sample sizes, may be read off from charts prepared 
by Clopper and Pearson (1934). Also Wilks (1943) has shown that 
approximate confidence limits of p, , corresponding to the confidence 
level a/100 may be obtained by solving the following quadratic in p, : 


where u is the normal deviate corresponding to the a% level of signifi- 
cance. On the other hand if p, is being estimated by the system of 
sampling outlined in section 4, its approximate confidence limits may 
be obtained as follows: 


Consider the function d = m, — p.m. , 
we have E(d) = 0; 


and V(d) = — [V(m.;) — 2p, Cov(m,; , + p:V(n,.)]. (18) 


If K is effectively large, then the factor (K — k)/(K — 1) may be 
taken as unity. In large scale sampling, m, and n, will be distributed 
almost normally. Hence 


u = — pan) [V(m,,) — 2p, Cov(m,; , ni) + (14) 


will be distributed almost normally with zero mean and unit variance. 
Hence the confidence limits of p, , corresponding to any desired. confi- 
dence level, may be obtained by solving the quadratic (14) for p, . 
However, since the population values of the covariance and the variance 
terms in this equation are unknown, the sample values are the only 
one available for use. Hence the sample size should be large to pre- 
serve the normality of u, and preferably k > 150. The principle of this 
method is due to Fieller (1940). This particular problem is just one 
aspect of the general problem discussed recently by Finney (1949) in 
this Journal. It is presented here in the notation used in section 4 to 
show that it has applications in demography. 
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THE USE OF SYSTEMATIC 5 X 5 SQUARES’ 


@IVIND NISSEN? 


paeeaang SQUARES were once used quite commonly by research 
workers; one of these, the 5 X 5 square which Fisher’ (3, page 76) 
calls the Knut Vik square, is still frequently used in Scandinavia. There 
are two possible different arrangements of this design, one with treat- 
ment A below C, the other with A below D. 


I II 
A BC ODE 48 ¢ 
DEABC c 8 
B C DEA EAs DBD 
E A BC D BCODEA 
C D E A B DE ABC 


The statistital analysis of these designs were originally made by 
adjusting for row and column differences. Since Fisher’s analysis of 
variance technique became known, this technique, which in effect gave 
the same result as the old method, was commonly used. From the 
total SS the SS due to rows, columns, and treatments were deducted, 
and the rest was regarded as error SS. Fisher (l.c.) showed that this 
method of analysis gave a biased estimate of the error. With a Latin 
square the calculated error is an unbiased estimate of the true error. 
If, by use of the proper systematic square, however, a lower true error 
is obtained, the calculated error will automatically increase. Tedin (4) 
confirmed this statement by calculations on uniformity trials. He 
found that the true error variance by use of a Knut Vik square was on 
the average 91.2 per cent of the error variance for Latin squares. The 


1Contribution from the Department of Agronomy, The Pennsylvania State College, State College, 
Pa. Authorized for publication on December 3, 1948, as Paper number 1496 in the Journal Series of 
the Pennsylvania Agricultural Experiment Station. 

2Visiting Professor, Department of Agronomy, Pennsylvania State College, on leave from Farm 
Crop Institute, Agricultural College of Norway, Vollebekk, Norway. 

3Numbers in parentheses indicate Literature cited. 
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calculated error variance was found to be 13 per cent greater than the 
true error variance. 

Dorph-Petersen (2), however, has shown that it is possible to calcu- 
late an unbiased estimate for the error in a Knut Vik square. For an 
ordinary Latin square there are two restrictions; each treatment shall 
appear once in each row and once in each column. Now, if two new 
restrictions are added, namely that each treatment shall appear once 
in each diagonal (both ways) then one of the two possible Knut Vik 
squares will be obtained. All plots with the same number in the fol- 
lowing scheme are then to be regarded as belonging to the same diagonal: 


1 2 3 4 5 
5 1 2 3 4 
4 5 1 2 3 
3 4 5 1 2 
2 3 4 5 1 


Since rows, columns and both diagonals are orthogonal to each other 
and to the treatments, the correct analysis of variance for a Knut Vik 
square will have this form: 


Degrees of freedom 
Rows 


4 
Diagonals 1st way 4 
Diagonals 2nd way ........ 4 
4 
4 

24 


Error 


However, it is not necessary to carry out the whole analysis in order 
to find the error. With arrangement J (p. 1) it is necessary only to 
calculate the sums of the plots with the same letter in arrangement JJ. 
The variance between these sums then becomes the error variance. 
This procedure utilizes the fact that the design is a completely orthogonal 
5 X 5 square (Fisher, l.c. p. 81). This method of calculating the error 
in a systematic 5 X 5 square was proposed by Vik (5) as early as 1924, 
and has been called “Vik’s short-cut method”. Dorph-Petersen (l.c.) 
confirmed these theoretical considerations by comparing F-values from 
uniformity trials with the theoretical F-distribution and obtained a 
satisfying agreement, at least for P = 0.05. 

The main disadvantage of the systematic square is that there are 
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only 4 degrees of freedom for error whereas there are 12 ina 5 X 5 
Latin square. 

According to Cochran & Cox (1, p. 29) the relative efficiency for 
two designs with the same error variance and with n, and n, degrees 
_of freedom for error is 


(mn, + 1)(m2 + 3) 
(n2 + 1)(m + 3) 


When comparing a Knut Vik square (n, = 4) with a Latin square 
(n. = 12) this fraction is (5 X 15)/(13 X 7) or 82.4%. As mentioned, 
Tedin (l.c.) found that the average true error variance for the Knut 
Vik square was 91.2% of the error variance for the Latin square. For 
a single experiment the ordinary Latin square is therefore the most 
efficient. In many cases, there is no special interest in the error of a 
single experiment. When a series of experiments is conducted in different 
years and in different localities, it is necessary to determine an average 
value for the experimental error in order to test the significance of the 
different interactions (treatment X year, treatment X locality, etc.). 
In this case, the relative efficiency of the systematic square is much 
higher than for a single experiment. For a small series of 5 experiments 
nm, = 20 (Knut Vik square) and n, = 60 (Latin square) and the fraction 
(21 X 63)/(61 X 23) = 94.38%. If it is assumed that the error variance 
for the systematic square is 91.2% of the error variance for the Latin 
square, the systematic square will be the more efficient one. For larger 
series this advantage of the systematic square will be still more pro- 
nounced. 
In such cases the important error is the error “between experi- 
ments’. This error will be the composite of the true variability be- 
tween experiments and the variability within each experiment. The 
first part is independent of the design used. If therefore, by using the 
Knut Vik square instead of the Latin square the later part can be 
reduced, the error “between experiments” will also be reduced thus in- 
creasing -the precision of the whole series of experiments. This in- 
creased precision is obtained without any increased cost. In fact, by 
using the Knut Vik square, the cost of planning, handling, and super- 
vising a series of similar experiments will be lower than by using Latin 
squares, and the chance of making errors will decrease. By means of 
Vik’s short-cut method, the labor involved in the calculation of the 
error within experiments will also be less than by use of Latin squares. 
The only disadvantage is that to some extent it is necessary to sacrifice 
some information from individual experiments. In working with a large 
series of experiments, the information from a single experiment is, how- 
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ever, of minor importance and in such cases the use of the Knut Vik 
square can be recommended. Of course the analysis of data obtained 
with these squares depends on the assumption of normality; there being 
only two possible Knut Vik squares, one cannot use the randomization 
argument to justify F tests as Pitman and Welch have done for ordinary . 
Latin squares. On the other hand, results of Dorph-Petersen in com- 
paring F-values from uniformity trials to the theoretical F distribution 
provide justification for using F tests. 
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THE EFFECT OF SULFAMERAZINE ON THE ERYTHROCYTE 
AND HEMOGLOBIN CONTENT OF TROUT BLOOD 


JAMES 8. GUTSELL 


U.S. Department of the Interior, Microbiological Laboratory 
Fish and Wildlife Service 
Kearneysville, W. Va. 


INTRODUCTION 


i= DISCOVERY, at Leetown in 1945, of the effectiveness of sulfa- 
merazine in the treatment of furunculosis in trout (1, 2), and the 
consequent widespread use of this drug for this purpose, made it im- 
portant to learn the effects of this sulfonamide on trout. Accordingly 
an experiment was run to test the effects of various dosage rates, from 
approximately the minimum effective rate to one greatly in excess of 
the highest rate recommended. 

The-daily dosage rate originally recommended (2) was 8 grams per 
100 pounds of fish. Later it was found that a lighter rate was sufficient 
and that a 4-gram rate was about as effective as the 6-gram rate subse- 
quently recommended (3, 4, 5). The dosage rates used in this experi- 
ment were 0, 5, 10 and 15 grams of sulfamerazine per day per 100 
pounds of fish. 

The duration of treatment originally recommended was three weeks. 
In this experiment sulfamerazine was administered through more than 
twice this period. 

In the original experiment (1, 2) the species of trout used was the 
Eastern brook trout, Salvelinus fontinalis, a species very susceptible to 
furunculosis. This experiment, described in part by Gutsell and Snieszko 
(12), included two other species: the brown trout, Salmo trutta or Salmo 
fario, generally less susceptible; and the rainbow trout, Salmo gairdnerii, 
which only occasionally suffers noticeable loss from this disease. 

The test using brook and rainbow trout was run in 1946; that using 
brown trout in 1947. Erythrocyte counts and hemoglobin measure- 
ments were made only in 1947 and therefore could include only the 
brown trout. 
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With rainbow trout, growth was not affected even by the heaviest 
dosage. With brook trout, growth was much reduced only by the 
heavier dosage rates; with brown trout, even by the lightest rate. . The 
reduced growth in brook and brown trout appeared to be due to reduced 
food consumption due to dislike of the taste of the food containing the 
drug. In this and in other experiments run at Leetown, there was not 
any mortality attributable to the use of sulfamerazine. 


METHODS AND PROCEDURE 


As noted above, the daily sulfamerazine dosage rates were 0, 5, 10 
and 15 grams per 100 pounds of fish (0, 11, 22, and 33 grams per 100 
kilograms). The drug was administered by being mixed very thoroughly 
with the food. Except for the 0-rate, all food used in the experiment 
contained the drug. Treatment was continued through 49 days. 

Five pounds of brown trout from one large pool were placed in each 
of 8 adjacent upper troughs in the hatchery. Each of the 4 treatments 
was assigned to 2 troughs by random selection, so that the test was 
run in duplicate and was fully randomized. 

The hemoglobin measurements were made after 35 days of treat- 
ment; the erythrocyte counts, after 42 days. A Helige hemometer was 
used for the hemoglobin measurements and a Spencer Bright-line count- 
ing chamber for erythrocyte counts. The usual technique for these 
operations was followed, except that Dr. S. F. Snieszko, who made the 
measurements and counts, used his cold method (Pipettes kept at a low 
temperature and fish placed in ice-water, to prevent the clotting which 
is the bane of such work with trout blood). For the erythrocyte counts, 
a 1-100 dilution was used. 

The fish used for the hemoglobin measurements and red cell counts 
were taken at random from each trough by means of a dip net. 

For hemoglobin measurements, 10 fish were taken from each trough, 
or 20 fish for each treatment. For erythrocyte counts, only half this 
number of fish were used, but two counting chambers (a and b) were 
filled from each fish. Thus there were 20 erythrocyte counts and 20 
hemoglobin measurements for each treatment. 


THE DATA AND THEIR ANALYSIS 


The average hemoglobin and erythrocyte contents are shown in 
Figure 1; the hemoglobin measurements and erythrocyte counts in 
Tables 1 and 3. Hemoglobin is expressed as grams per 100 milliliters 
of blood. Each erythrocyte count in Table 3 is the sum of the counts 
for 5 secondary squares. Because the dilution was 1-100, each count 
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in Table 3 (or each average of such counts) should be divided by 2 
before the 4 zeros (0, 000) are added to give the number of erythrocytes 
per cubic millimeter of blood.’ 
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1 1 1 1 
5 10 15 
DOSAGE 
FIGURE 1. 


Hemoglobin (grams per 100 ml.) and erythrocytes (millions per cu. mm.) 
in the blood of brown trout after 35 (hemoglobin) and 42 (erythrocytes) 
days of treatment at the daly dosage rates of 0, 5, 10 and 15 grams of 
sulfamerazine per 100 pounds of fish. Each point represents the average 
of 20 determinations. 


1Each count is the number of red cells in 5 secondary squares of the counting chamber. The 
volume of the diluted blood in the 5 secondary squares is 1/50 of a cubic millimeter. Therefore, the 
estimated number of red cells in a cubic millimeter of the diluted blood is 50 times the count. If the 
dilution was 1/100, the estimated number of erythrocytes per cubic millimeter of blood is 5,000 
(100 X 50) times the count. This multiplication is most conveniently done by dividing the count by 
2 and adding 4 zeros (0,000; equivalent to multiplying by 10,000). 
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TABLE 1 
MEASUREMENTS OF HEMOGLOBIN (GRAMS PER 100 ML.) IN THE BLOOD OF BROWN 
TROUT AFTER 35 DAYS OF TREATMENT AT THE DAILY RATES OF 0, 5, 10 AND 15 
GRAMS OF SULFAMERAZINE PER 100 POUNDS OF FISH. TWO TROUGHS, A AND B, 
FOR EACH TREATMENT; 1 MEASUREMENT FROM EACH OF 10 FISH FROM EACH 


TROUGH. 
Treatment 0 Treatment 5 | Treatment 10 | Treatment 15 
Fish no. 
A B A B A B A B 
rf 6.7 7.0 9.9 9.9 10.4 9.9 9.3 | 11.0 
2 7.8 7.8 8.4 9.6 8.1 9.6 9.3 9.3 
3 5.5 6.8 | 10.4 10.2 | 10.6 10.2 7.2 11.0 
4 8.4 7.0 9.3 10.4 8.7 10.4 7.8 9.0 
5 7.0 7.5 10.7 | 11.3 10.7 | 11.3 9.3 8.4 
6 7.8 6.5 11.9 9.1 9.1 10.9 10.2 8.4 
7 8.6 5.8 7.1 9.0 8.8 8.0 8.7 6.8 
8 7.4 7A 6.4 10.6 8.1 10.2 8.6 2 
9 5.8 6.5 B36. 141.7 7.8 6.1 9.3 8.1 
10 7.0 5.5 10.6 9.6 8.0 | 10.7 7.2 111.0 
Trough sums | 72.0 | 67.5 | 93.3 |101.4 | 90.3 | 97.3 | 86.9 | 90.2 
Treatment 
sums 139.5 194.7 187.6 TA 
means 6.98 9.74 9.38 8.86 
Total 698.9 


Analyses of variance of the data are given in Tables 2 (hemoglobin) 
and 4 (erythrocytes). The higher hemoglobin and erythrocyte contents 
of the blood of fish given sulfamerazine are found to be highly significant 
(F for treatments above the 0.1 percent probability points). For the 
erythrocyte data, this applies not only to the counts but also to the 
square root and logarithmic transformations (see bottom of Table 4). 
Neither transformation changed the picture. 

In the breakdown of the three degrees of freedom for treatment, for 
each degree F is above the 5 percent probability point, and for all but 
one degree (third degree, hemoglobin) above the ene precent probability 
point. 

To get at this another way, the values for the 5-gram rate were 
compared with those for each of the others. For hemoglobin the com- 
parisons were as follows: 5-gram rate vs. 0-rate, F = 49.11; 5-gram rate 
vs. 10-gram rate, F = 0.81; 5-gram rate vs. 15-gram rate, F = 4.99 
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TABLE 2 
ANALYSIS OF VARIANCE OF HEMOGLOBIN DATA. 
Degrees | Sums Probability Points 
Source of of of Mean F 
variation freedom | squares | squares 5% 1% 0.1% 
Total 79 209.48 
Treatments 3 90.56 | 30.19 | 19.48 2.74 4.08 6.09 
(1st°) Linear . 1 27.93 | 27.93 | 18.02 3.98. 7.01 11.87 
(2nd°) Quad- 
vratic ... 1 53.96 | 53.96 | 34.81 3.98 7.01 11.87 
(8rd°) Cubic . 1 8.67 | 8.67 5.59 3.98 7.01 
Between 
troughs : 
treated alike 4 7.29 1.82 2.50 
Between fish in 
same trough 72 111.64 | 1.55 
TABLE 3 


ERYTHROCYTE COUNTS! (5 SECONDARY SQUARES; DILUTION, 1-100) FROM BROWN 
TROUT AFTER 42 DAYS OF TREATMENT AT DAILY RATES OF 0, 5, 10 AND 15 GRAMS 
OF SULFAMERAZINE PER 100 POUNDS OF FISH. TWO TROUGHS, A AND B, FOR EACH 
TREATMENT; 5 FISH FROM EACH TROUGH; 2 SLIDES, a AND b, FROM EACH FISH. 
TREATMENT MEANS SHOWN (IN PARENTHESES) AS NUMBERS PER CU. MM. 


Treatment 0 Treatment 5 Treatment 10 Treatment 15 
Fish 
no. A B A B A B A B 
a b a b a b a b a b a b a b a b 
1 213] 230] 166) 157] 296) 319] 310| 309} 339] 322) 196) 232) 278] 212) 287| 280 
2 253) 231) 206) 185} 278] 258] 241; 270} 282) 285) 205) 186} 275) 311] 221) 243 
3 195) 164) 245) 250) 345) 307) 272) 311) 236) 262) 186) 158) 331] 309 
4 203} 213) 181} 322) 372) 254) 237] 252) 209) 245) 216} 301) 281] 244 
5 191} 195} 198} 169] 248) 266) 275} 263] 296) 249} 260} 223) 246) 292) 295 
Column sums 1045} 1023/1028} 942) 1489) 1530) 1343) 1402) 1372|1374| 1147/1168] 1263) 1208) 1362|1371 
Trough sums. . 2068 1970 3019 2745 2746 2315 2471 2733 
Treatment sums 4038 5764 5061 5204 
Treatment means 201.9 288 .2 253 .05 260 .2 
(1,009,500) (1,441,000) (1,265,250) (1,301,000) 
20,067 


1To convert counts to numbers per cubic millimeter, divide by 2 and add 0,000; thus 230/2 = 115; 
115 plus 0,000 = 1,150,000. 


(See footnote to text). 
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(F .05 = 3.98). For erythrocytes: 5-gram rate vs. 0-rate, F = 30.45; 
5-gram rate vs. 10-gram rate, F = 5.05; 5-gram rate vs. 15-gram rate, 
F = 3.21 (F .1 = 2.88, F .05 = 4.15). 

In a ¢ test, a comparison of the 5-gram rate with the 15-gram rate 
(hemoglobin) gave ¢ = 2.078 (¢ .05 = 2.025). 


TABLE 4 
ANALYSIS OF VARIANCE OF ERYTHROCYTE COUNTS 
Source of Degrees of | Sums of Mean Probability points 
Variation freedom squares squares Fi 
5% 1% 0.1% 
79 186,948 .89 
Treatments. ..... 3 77,667 .74| 25,889.25 | 10.58 2.90 4.46 6.96 
(ist?) . . 1 19,530 19,530 .06 7.98 4.15 7.50 13.15 
(2nd°) Quadratic 1 31,323 .61| 31,323.61 | 12.81 4.15 7.50 13.15 
(3rd°) Cubic ..... 1 26,814 .06| 26,814.06 | 10.96 4.15 7.50 13.15 
Between troughs treated 
4 16,954.25) 4,238.56 1.73 2.67 
Between fish in same 
32 78,267 .40| 2,445.86| 6.93 1.74 2.20 2.6 
Between samples from 
same fish. ..... 40 14,059 .50 351.49 


1After square root transformation, the values found for F for treatments, troughs and fish were 
11.77, 1.82, 4.37; after logarithmic transformation, 11.01, 1.61 and 4.01. 


DISCUSSION 


The analyses of the data strongly indicate that the higher hemoglobin 
and erythrocyte concentrations in the blood of the brown trout on the 
diets including sulfamerazine is causally connected with the presence of 
sulfamerazine in the diet. There is considerable indication that the 
higher content associated with the 5-gram drug rate is significant (see 
Tables 2 and 4, individual degrees of freedom for treatment, also 
reference in the text to tests of significance), that the 5-gram rate truly 
is superior in this respect (perhaps near the optimum rate) under the 
conditions of the experiment. 

An important consideration is the relationship of the hemoglobin 
and red cell concentrations associated with the various treatments, 
including the 0-rate, to the normal concentrations. What are the normal 
concentrations? Schlicher (6, quoted by Wunder, 7) listed 1,140,000 
erythrocytes per cubic millimeter of blood for brown trout, sensibly the 
same as for rainbow trout (1,100,000). Phillips (8, Fig. 1) found above 
1,500,000 erythrocytes per cubic millimeter of blood from brook trout 
before they were put on a synthetic diet and after they were given beef 
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liver to the extent of half of the diet. Tunison et al, (9) present these 
same findings and others showing brook trout on a straight beef liver 
diet with 1,176,666 red cells per cubic millimeter. Tunison et al, (10) 
reported the following yearly averages for trout on a diet of beef liver 
(50 percent), pork spleen (48 percent) and salt (2 percent): rainbow, 
1,220,000; brook, 1,032,000; brown, 766,000 erythrocytes per cubic 
millimeter of blood. Phillips et al, (11) found the following number of 
red cells per cubic millimeter of blood after the trout had been 22 weeks 
on a diet of beef liver (98 percent) and salt (2 percent): rainbow, 
1,165,000; brook, 1,186,000; brown 1,099,000. On a diet of beef spleen 
(98 percent) and salt, the numbers were very similar: rainbow, 1,164,000; 
brook, 1,181,000; and brown, 1,116,000. Schlicher (see above) gave the 
hemoglobin level for brown trout as 60 percent (presumably of the 
normal for human blood or .6(14.5) = 8.7 grams per 100 milliliters). 

If it were not for the one high level, shown by Phillips (8, Fig. 1) 
and Tunison et al, (9, Chart 11), a fair assumption from the above 
references would be that for brook, rainbow and brown trout the normal 
is close to 1,000,000—1,200,000 per cubic millimeter. This is above the 
average for the “control” trout in this experiment, below the average 
for the others. If about 1,500,000 (8, 9) is normal, only the trout 
treated at the 5-gram rate showed an average approaching the normal 
red cell concentration. 

Another important consideration is the probable sulfamerazine con- 
centrations or levels obtained with the various treatment rates. For 
brown trout, data are available for approximately the 5-gram rate. 
Snieszko and Friddle (unpublished) used the daily rate of 0.1 gram of 
sulfamerazine per kilogram of trout (4.54 grams per 100 pounds) with 
this species. After a week of treatment the following levels were found: 
blood, 64; liver, 24; kidney, 24; muscle, 11 milligrams percent. There- 
after, the following more moderate levels (approximate) prevailed: 
blood, liver and kidney, 8; muscle 5. This was the only rate studied 
with brown trout; but with brook trout, 0.1, 0.2 and 0.3 grams per day 
per kilogram of trout (roughly corresponding with the 5, 10, and 15 
grams per 100 pounds used in this experiment) were investigated. 
The following levels approximate those which prevailed after three days 
in gutted fish and presumably are close to the levels for the muscle: 
0.1, 3; 0.2, 6; and 0.3, 7 milligrams percent. These findings indicate that 
the 5-gram rate is unusually effective with brown trout in producing 
satisfactory therapeutic sulfamerazine levels, and that the high hemo- 
globin and erythrocyte levels which followed prolonged treatment with 
sulfamerazine were associated with fairly high sulfamerazine levels in 
the blood and tissues. 
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That the sulfamerazine treatment, almost certainly accompanied by 
high sulfonamide blood levels, did not reduce the hemoglobin and 
erythrocyte concentrations is notably important; that these concentra- 
tions were increased is of great interest and of possible importance. 
The considerably higher concentration associated with the use of one 
of the three dosage rates is of additional interest. 

Because prolonged sulfamerazine treatments at high dosage rates 
had not been accompanied or followed by any trout mortality attrib- 
utable to the drug (2, 5, 12), the lack of a reduction in the concentration 
of hemoglobin and erythrocytes was not altogether surprising; but the 
increase was a great surprise and still is a puzzle. Gutsell and Snieszko 
(12) offered the following possible explanations: 1, some infection, ap- 
parently not furunculosis (13), causing a moderate anemia was cured 
or greatly reduced by sulfamerazine therapy; 2, the sulfonamide to 
some degree inhibited the respiratory enzymes (14) with the effect of a 
reduction in the oxygen supply which is (15) to increase the hemoglobin 
and erythrocyte concentrations; 3, sulfamerazine, instead of having an 
anemic effect, stimulates the production of hemoglobin and erythrocytes 
in trout. The writer offers a fourth possible explanation, that the diet 
caused a moderate anemia, and that the sulfamerazine, by reducing the 
consumption of food (as was observed to be the case), reduced the 
anemia (and growth, of course). 

The clear evidence that continued dosage with sulfamerazine did not 
reduce the hemoglobin or erythrocyte concentrations explains in some 
measure why adverse effects, other than poor food consumption (ap- 
parently due to dislike of the taste of the drug) and reduced growth, 
have not been observed at Leetown, either in experiments or in practical 
application of sulfamerazine therapy for furunculosis. 

The author wishes to thank Dr. C. M. Mottley for suggestions for 
the analysis of these data, especially the erythrocyte counts, and for 
reading the manuscript. Dr. S. F. Snieszko, as a collaborator in the 
investigation, made the erythrocyte counts and hemoglobin measure- 
ments. Mr. Saufley Friddle helped in the preparation of the diets and 
in feeding the trout. Messers. E. W. Surber and B. H. Hazen made 
hatchery troughs available and extended numerous courtesies. 
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SUR UNE APPLICATION DE LA STATISTIQUE 
MATHEMATIQUE A LA BIOLOGIE 


Maurice Fr&cHET 


a Paris 


Rappel.—Nous avons consacré ces derniéres années trois mémoires' 
4a l’étude des éléments aléatoires de nature queleonque, en vue de 
généraliser la théorie des nombres aléatoires, objet essentiel, jusqu’d 
présent, du calcul des Probabilités. Dans cette étude, nous avons 
défini de plusieurs maniéres les positions typiques d’un élément aléatoire 
de nature quelconque, par extension de la notion de “valeur typique”’ 
(valeur moyenne, équiprobable, dominante, etc.) d’un nombre aléatoire. 
Si, par exemple, on a défini une “‘distance’’ (a, b) de deux quelconques, 
a, b, des déterminations que peut prendre |’élément aléatoire X, on 
pourra opérer ainsi: Soit @(a) la racine d’ordre a de la valeur moyenne 
de la puissance a(>0) de la distance (X, a) 


= VMN(X, a)’. 
8’il existe une position y, de |’élément certain a, tel que l’écart moyen, 
6(a), d’ordre a de X avec a, atteigne sa borne inférieure quand a 
devient y. , nous appellerons y, une position typique d’ordre a de X. 
' Ce sera une position moyenne pour a = 2, une position équiprobable 
pour a = 
Quand, par exemple, |’élément X est une fonction X(é) choisie au 


hasard parmi une famille F de fonctions numériques continues f(f) sur 
(0, 1) et si 


(f, 9) = max | f(t) — g(t) | 


la position moyenne de cet élément, X(b), est une fonction continue 
qui n’a pas une définition particulitre pour chaque valeur de ¢ mais 


1M. Fréchet, Les éléments aléatoires de nature queleonque dans un espace distane:é, Annales de 
UInstitut H. Poincaré, 14 (1948), p. 215-310; L’intégrale abstraite d'une fonction abstraite d'une 
variable abstraite et ses applications 4 la moyenne d’un élément aléatoire de nature quelconque, 
Revue Scientifique, 82 (1944), p. 483-485; Positions typiques d’un élément aléatoire de nature quel- 
conque, Annales Ecole Normale Supérieure, 65 (1948), p. 211-237. 
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qui est définie comme un bloc. Elle doit étre distinguée de la fonction 
X(é) qui, pour chaque valeur de ¢, est la moyenne au sens classique du 
nombre aléatoire X(#). Et nous avons montré par un exemple que ces 
deux définitions peuvent conduire 4 des fonctions numériquement 
différentes. 

D’une facgon générale, quand un élément aléatoire Y est déterminé 

par la connaissance de certaines caractéristiques numériques Y, , Y,, - 
il n’y a pas de raison pour que notre définition conduise & prendre pour 
élément moyen Y un élément déterminé par des caractéristiques y, = 
Y:, ¥ = Ys, +--+ respectivement égales aux moyennes de Y, , de 

Application Vhomme moyen.—Ayant eu tout récemment l’oceasion 
d’entendre renouveler les critiques—d’ailleurs fondées—adressées 
maintes fois 4 la conception, due 4 QuETELET, de l’homme moyen, il 
nous a paru qu’on pouvait les éviter en modifiant la définition de l'homme 
moyen conformément 4 la derniére remarque ci-dessus. 

En effet, d’aprés QureTELET, l’-homme moyen (d’une population 
donnée) serait par définition celui dont chaque caractéristique (poids, 
taille, etc.) serait mesurée par la moyenne des valeurs de cette carac- 
téristique (dans la population donnée). 

Joseph BERTRAND a fait observer qu’une telle définition pourrait 
conduire 4 un résultat impossible ou anormal, que, par exemple, le 
poids moyen pourrait n’étre le poids d’aucun individu ayant la taille 
moyenne ou (plus vraisemblablement) étre beaucoup plus fort que les 
poids de la plupart des individus ayant la taille moyenne. 

On échappera 4 cette objection de principe en appliquant 4 une 
population humaine notre définition générale ci-dessus. Nous l’expo- 
serons sous une forme un peu modifiée. Supposons qu’on puisse, con- 
naissant deux individus J, J’ et leurs caractéristiques p, t, «++ , p’, l’, --: 
estimer numériquement leur dissemblance au moyen d’un nombre 
(I, I’) qui sera une fonction déterminée des différences en valeurs 


absolues |p — p’|, |¢ — ¢’|. - - - (On pourra, par exemple, prendre 
pour (J, I’) une valeur typique de l’ensemble des nombres | p — p’ |, 
|¢ — t’|, ---). Alors, on pourra définir un degré numérique de dis- 


semblance de J avec |’ensemble de la population P au moyen d’une 
valeur typique [J, P] de l’ensemble des valeurs de (J, I’) pour I fixe, I’ 
variable pour P. 

Ceci étant, l’-homme typique (de la popuiation) sera celui, J, dont 
le degré de dissemblance [J, P| avec la population sera le plus faible, 
c’est-a-dire tel que [J, P] < [J, P] quel que soit l’individu J de la 
population P. 

Remarque: I arrivera généralement que l’homme moyen (ou typique) 
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selon la nouvelle définition sera exceptionnel par quelque trait. II 
échappera ainsi 4 la critique de BerTRAND selon laquelle l’homme 
moyen serait “‘médiocre en tout”’. 

Si l’on enléve 4 la critique de BERTRAND sa pointe (car ici “‘médiocre”’ 
n’est que le sens péjoratif de ‘‘moyen’’), on peut cependant se demander 
si elle concerne vraiment un défaut et s’il ne conviendrait pas que 
homme moyen soit,—sinon moyen en tout comme le voulait QUETELET 
ce qui peut étre impossible—, du moins qu’il soit aussi prés que possible 
d’étre moyen en tout. C’est ce qui nous a amené, indépendamment, 
Mme Derrise (4 qui j’avais communiqué la définition ci-dessus) et 
moi-méme a définir autrement l’homme moyen: ce serait ainsi celui des 
membres de la population considérée dont l’ensemble des caractéristiques 
différe le moins de l’ensemble des caractéristiques moyennes. La 
différence globale pourra étre mesurée par une valeur typique de 
l’ensemble des valeurs absolues des différences entre chaque caractér- 
istique de l’individu considéré et la caractéristique moyenne corre- 
spondante. On pourra appeler l’homme ainsi choisi, l’homme le moins 
éloigné de la moyenne ou plus briévement l’homme le plus moyen. On 
définirait, de méme, l’homme le plus médian, --- l’homme le plus 
typique. 

Nous développerons ailleurs plus en détail ces diverses définitions 
de l’homme typique par exemple en ce qui concerne le repérage numé- 
rique des caractéristiques. 

Détermination de l’espéce-——Nous avions remis en Belgique puis 
envoyé de Paris deux notes successives sur cette conception de l’homme 
moyen ou typique (pure application de nos conceptions générales anté- 
rieures, rappelées ci-dessus), lorsque nous avons lu & notre retour 4 
Paris, un article, court mais trés intéressant, par M. S. KuiczyXsk1 
sous le titre “Sur quelques applications des notions topologiques a la 
biologie’’.” On verra tout de suite la relation de la proposition de cet 
auteur avec la nétre, si on la résume en deux points: 1°. On représente 
chaque individu J d’une population animale ou végétale par un certain 
nombre de valeurs caractéristiques P, 7, --- qui peuvent étre con- 
sidérés comme ses coordonnées dans dans un certain espace. L’ensemble 
des points correspondant aux individus de la population forme un 
nuage; 2° Chaque “point d’accumulation” de ce nuage détermine un 
individu de la population qui définira une espéce. 

Ainsi nous voyons que chaque espéce est représentée par un individu 
pris dans l’espéce elle-méme, ce qui le fait échapper 4 l’objection de 
BERTRAND, et qui constitue ainsi une amélioration de la définition de 


2Colloquium Math., 1 (1948), p. 182-183. 
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QUETELET. Bien que l’auteur considére qu’il y a 1a une application de 
la topologie, nous croyons qu’il s’agit plutdt d’une application de la 
statistique mathématique. L’auteur laissant indéterminée la définition 
de ce qu’il appelle un point d’accumulation, nous pensons que son inter- 
prétation la plus satisfaisante serait une position ot la densité statistique 
atteint un maximum local’ ou mieux comme une position typique au 
sens de la statistique mathématique sous la forme généralisée ci-dessus 
pour l’un des nauges de points séparés par des régions 4 densité statis- 
tique nulle on minime. 

Généralisations.—I] nous a paru en tout cas intéressant de signaler 
la convergence (en méme temps que l’indépendance) des conceptions 
de M. KuuczyXskr et des nédtres. Et aussi de faire observer que ces 
conceptions de |’individu représentatif d’une espéce ou d’une population 
d’étres vivants peuvent étre étendues au cas beaucoup plus général 
d’une population d’éléments de nature quelconque. C’est ainsi que 
nous rassemblons actuellement les éléments statistiques permettant de 
définir une ville typique (d’un ensemble assez homogéne de villes 
d’importances comparables). 

Dans les mémoires cités plus haut®* nous avons d’ailleurs aussi 
étudié d’autres secteurs (comme, par exemple, |’étude des formes 
variées de convergence stochastique) de la généralisation de la théorie 
des nombres aléatoires 4 celle des éléments aléatoires de nature quel- 
conque. 

Addendum.—Lors de notre exposé sur ce qui précéde, devant la 
nouvelle Société frangaise de Biométrie, le 17 mai 1949, un de nos audi- 
teurs, M. DreLaporte, nous a signalé qu’il avait, indépendamment de 
M. KutuczyNsk1, exposé des idées analogues 4 celles de ce dernier, en 
1946, aux pages 14-23 du tome 8 du périodique édité par la Société de 
Biotypologie (et méme, auparavant, dans une communication non 
publiée faite au Congrés de l Association francaise pour 1 Avancement des 
Sciences en Octobre 1945). On trouvera une mise au point récente de 
ses idées dans son article “Bases théoriques de la détermination des 
types’’, dans la Revue Sciences, Paris, 76(1949), p. 110-119. 

Nous nous étions rencontré avec lui dans notre communication en 
remplacant le point d’accumulation de M. KunczyNski déterminant 
l’espéce (ou le type) par un point de densité maximum locale qu’on 
peut appeler point dominant local. 

Bien entendu, il n’y a pas lieu, généralement, de considérer ‘tout 
point ot la densité de probabilité (ou, en pratique, la densité de 


3(Cf. note !). Nous avons résumé une partie de ces mémoires: en francais dans le Journal de la 
Société de Statistique de Paris, 88 (1947), p. 410-420 et en anglais dans Mathematics Magazine, 22 
(1948-49), p. 1-12. 
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fréquence) est plus grande que pour les points voisins, comme définissant 
un type ou une espéce. II pourra arriver que ce soit simplement un 
point remarquable de la répartition de la population relativement & ses 
caractéres ou simplement un effet du hasard. II faudrait s’assurer si l’on 
peut répartir la population en plusieurs lots dont les nuages de points 
représentatifs correspondants s’amassent nettement chacun autour 
d’une position déterminée. Et celle-ci ne sera pas nécessairement un 
point dominant mais un point moyen ou médian, --- ou en général 
typique au sens précisé plus haut. La distinction en types sera plus 
nette si les nauges correspondants sont complétement séparés les uns 
des autres. Mais cela n’arrivera guére que pour des types trés différents 
ou dans les cas ot la population totale est peu nombreuse, cas ov 
l’existence de nuages séparés pourra étre due au hasard. D’ailleurs la 
notion de nauges séparés n’est pas aussi claire qu’il parait d’abord. En 
résumé, méme au seul point de vue théorique, la notion de type ou 
d’espéce est devenue, gréis 4 ce qui précéde, plus précise, sans étre 
encore complétement élucidée. 
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RELATIONSHIPS AMONG CHARACTERS OBSERVED IN 
UREA CLEARANCE TESTS 


H. FarrFIELD SMITH 


Institute of Statistics 
North Carolina State College 


I LABORATORIES of clinical pathology the observation usually made 
to measure efficiency of kidney action is that known as the van Slyke 
urea clearance test. Failure of renal function results eventually in 
increased concentration of urea in the blood; the objective of the test 
is to provide a measure of partial failure before the onset of that conse- 
quence. 

It has been supposed that this can be obtained by observing the 
rate of urea excretion in urine. Following observations (under some- 
what limited and special conditions) that urea excretion is approxi- 
mately proportional to concentration in the blood, the measure adopted 
for test purposes has been the ratio of urea excreted per minute/urea 
concentration in blood—usually described as measuring the volume of 
blood “cleared” of urea per minute, although of course no blood is 
actually cleared. 

A second consideration which the van Slyke test was devised to meet 
is that rate of urea excretion depends also on the rate of water excretion, 
particularly when this is slow. It was observed that rate of urea excre- 
tion is maintained nearly constant when rate of urine excretion exceeds 
about 2 ml. per minute. Consequently “maximum urea clearance” was 
defined as 


C,, = UV/B (V > 2) (1) 


where V = volume of urine excreted in ml. per minute 

U = concentration of urea in urine measured as mgm. N per ml. 

B = concentration of urea in blood measured as mgm. N per ml. 
In their early work van Slyke and his colleagues observed that at lower 
rates of urine excretion the quantity of urine urea (UV) varied ap- 
proximately in proportion to~/V. Using this relationship to estimate 
quantity of urea excreted when V = 1, they defined “standard urea 
clearance’, 


C,=UVJV/B (V <2) (2) 
185 
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as an estimate of the ratio UV /B at a urine excretion rate of one ml. 
per minute. 

Text books of clinical methods usually quote formulae (1) and (2), 
with the accepted criteria for each, but with no indication of their 
derivation. This was the only information available when in a prisoner 
of war camp in Siam an attempt was made to run some such tests as 
described by Wilshaw, 1946. The men observed were not expected to 
be renally defective, having been sent for test only as a precaution prior 
to operation. At first sight the results appeared alarmingly low, but 
it soon became apparent that they could not be interpreted according 
to the usual criteria. Therefore after checking chemical procedure, a 
group of unquestionably normal men were observed with intent to 
develope our own standards. Having no clue to the reason for the 
two usual formulae it seemed advisable to derive also our own adjust- 
ment for rate of urine excretion. For this purpose a linear relationship 
between log U and log V was found satisfactory. 

When access to literature was again possible it seemed of interest to 
see if other published data might indicate that the method of adjust- 
ment used in Siam, or a similar one, could be generally adopted to 
provide a single criterion for any rate of urine excretion in place of 
rather artificial partition of observations into two groups with different 
criteria. Using the regression given below that can be done; but in 
the course of the work two fallacies underlying the test became apparent. 


1. The correlation of urine and blood ureas which provides the rationale for the 
test was observed only for response within individuals to temporary experi- 
mentally imposed alterations of blood urea concentration. Observations on 
prisoners of war on deficient protein data, corroborated by other observations 
in India, showed that as between individuals on permanently different diets 
this correlation is not operative. Consequently the comparison to blood urea 
does not make allowance for differences in diet, as would often be necessary for 
validity of the test as usually applied. 

2. The supposedly established criteria for normality were based on inadequate 
sampling procedure. They were merely means and ranges observed for a 
haphazardly selected group of 19 individuals (presumably workers around the 
laboratory). Such a sample gives little information on variation in the popu- 
lation at large. Furthermore, patients suspected of renal deficiency are 
commonly put on low protein diet and of such these nineteen are no sample 
at all. 


These aspects have been discussed in a separate paper (Smith, 1949). 
Since the basis of the test is thus open to serious question the original 
objective of the statistical work here reported is irrelevant for im- 
mediate application. However, in so far that a statistical summary of 
existing data may aid further investigation, it seems advisable to make 
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available a brief account of the work which has been done. Should it 
still be considered desirable, for any purpose, to make comparisons of 
urine urea concentrations (U) or urea excretions (UV), either for the 
same person at different times or between persons, a method is suggested 
by which observations may be adjusted for rate of urine excretion, for 
alterations of blood urea due to temporary fluctuations (but not for 
permanent differences), and for size of person (required only when 
persons of appreciably different size are to be compared). But it is 
suggested that such comparisons be not considered as tests for renal 
defect, at least until the effects of diet have been more thoroughly 
elucidated. 


INTERNAL REGRESSIONS 


The first objective was to establish the relation between urea con- 
centration and volume of urine. Figure 1 shows a typical set of data 
for one individual at approximately constant blood urea concentration. 
The procedure of van Slyke and his co-workers was equivalent to ap- 
proximating the relationship by two straight lines 


log (U/B) = log C, — 3 log V (V < 2) (3) 


log (U/B) = log C,, — log V (V > 2) (4) 


Their point of intersection has been called the “augmentation limit”, 
and has been discussed as a feature of potential interest to characterise 
different persons. As shown by Figure 1 these lines are not unreasonable 
as a first approximation, and (4) may be regarded as an asymptote 
for high V. It may be assumed, however, that the true curve should 
have no discontinuity; nor can (3) be regarded as asymptotic to the 
other extreme because all adults on whom more than ten observations 
have been made show that slope of log U on log V becomes numerically 
less than 0.5 when V is less than 1 ml. per minute (as has been noted 
also by Chesley, 1937, 1938, and by van Slyke, 1947). Consequently 
the “augmentation limit” is an artifact which cannot be regarded as a 
suitable statistic by which to describe curves exhibited by different 
persons. 

The true form of relationship being unknown’, available data can be 


*A curve asymptotic to (4) would appear to have some theoretical advantage. Several such were 
tried of which the best appeared to be (figure 1) 
log U = A — log (a + V) + k log B. (5) 
However, in the absence of knowledge of the true relationship, such forms appeared to give little 
advantage over (5), which is simpler for statistical purposes although unrealistic for large values of V 
beyond the point where du/dv = m + 2cv = —1. Within observable range such defect is trivial 
(figure 1). 
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FIG. 1—OBSERVATIONS ON JHA, DATA FROM AUSTIN ET AL (1921) 

Curve (3): log (U/B) = log 49.8 — } log V 
Curve (4): log (U/B) = log 77.2 — log V 
Curve (5): log (U/B) = 1.957 — log (.80 + V) 
Curve (6): log (U/B) = 1.710-.5191 log V — .2865 (log V)2 
(3) and (4) represent curves fitted by van Slyke et al, and correspond to Cs = 49.8, Cm = 77.2, respec- 
tively. This subject showed k =~ 1. 

Variances of deviations of log U from the regressions are: 
for (3) and (4) .00127, for (5) .00113, for (6) .00104. 


empirically described over the observable range and within limits of 
random variation by a polynomial of form. 


u=atm-+eo° + kb (6) 


where, for brevity, we write v = log V, u = log U, b = log B. Common 
logarithms have been used. 

All three variables, U, V and B being observed “‘as they come’’, in 
trying to evaluate the relationship between them there arises the ques- 
tion, as posed by Winsor (1946): ‘Which regression?” Throughout this 
work I take the point of view that our purpose is to examine what 
values of U occur with given observations of V and B (cf. postscript). 
Consequently I adopt the procedure of treating U as the dependent 
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variate and evaluate its regression on V and B as independent variates.** 
Although no definite evidence is available, it seems that random varia- 
tion may be roughly proportional to magnitude of obser vation; therefore 
log observations are taken to have equal weight. 

To seek evidence for the possibility that regression on blood urea 
might vary with rate of excretion (or conversely), the equation 


u=atm +o’ + kb + g(bv) 


was fitted to data for four men and for dogs reported by Austin ef al 
(1921). No evidence for that possibility was detectable. 

A separate regression of form (6) was fitted for each of the ten subjects 
for whom Austin et al (1921) and Moller eé al (1928) report more than ten 
observations. (Actually two sets of observations, treated here as if for 
separate subjects, were on the same person in two widely separated 
periods, 1921 and 1928.) Nearly all drank varying amounts of water to 
produce variation in rate of urine excretion. Five (1928) ate urea to 
increase blood urea content on some days. There was some indication 
that, apart from effects of urea consumption, blood urea concentration 
may have tended to decrease with advancing time of day and when much 
water was drunk; that it varied more from day to day than within days, 
and was nearly as variable for one person between days as between 
persons. It would be of interest to discriminate between normal varia- 
tion and superimposed effects of each factor; but as the experiments 
were not planned so as to facilitate such a breakdown it is scarcely prac- 
tical to attempt it. 

As a preliminary to further analysis we want to know whether or 
not all sets of observations show similar random variability. For 
variances of deviations about each regression, Bartlett’s test (with nota- 
tion of Thompson and Merrington, 1946) gives: 


for all ten variances M = 27.85, P < .01, 


excluding 2 subjects (Add and McL), k = 8, M = 12.22, 
P slightly > .05 


excluding 3 subjects (Add, McL and WN) k = 7, M = 6.18. - 


The data for Add and McL had been previously reported from other 


**It seems not to be usually recognized, and should therefore be pointed out, that even if we 
ascribe errors to V and B and fit a curve by the procedure described by Edwards Deming (1943), 
ignoring second and higher powers of residuals, we still arrive at the same curve as with the simpler 
assumptions. 
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laboratories and may therefore have been subject to different experi- 
mental variability. There is less justification for separating out WN, 
but some other features of his data also appear suspicious. If the 
statistical analysis were done on the spot one would like to repeat 
observations on this subject to see if he would persist in slightly ex- 
ceptional behaviour. In the following, data for Add have been weighted 
by 3, those for McL and WN by 3, these weights being, in round figures, 
approximately proportional to reciprocals of the ratios of the individual 
variances to the mean variance of the other seven. 


TABLE 1 
Variance 
Coefficient Range Weighted ratio, F. P 
mean* (m = 9, n. = 127) 

u’ (atv = .3, 

b = —.7) -70to .94) .815 + .022 7.00 <.001 
m’ (at v = .3) —.79 to —.60}—.701 + .019 2.28 ca .025 
c — .62 to —.16|—.355 + .050 3.87 <.001 
k .27to 1.14) .802 + .081 2.98 < .005 


*Standard errors determined from mean squares between persons. 


Figure 2 shows the fitted regressions on v, and table 1 summarises 
variation of the regression coefficients as computed for the form 


u=u' +m'(v — + — .3)? + k(b + .7) (7) 


That is, we compare urine urea concentrations (u’), and their rates of 
change with »v, (m’), near the middle of the range of observations, namely 
at v = .3, b = —.7. (Over-all averages are) = .313, b = —.755. 
If wu’ and m’ were evaluated for a non-central origin as in (6) they would 
be subject to higher errors. Coefficients c and k are not affected by 
change of origin for either variate, m’ is affected only for v.) The last 
two columns show tests of significance for heterogeneity of each type 
of coefficient after fitting 30 others (each of the other three for every 
individual). The mean squares for each coefficient are not independent; 
a test for heterogeneity of m, c, k collectively is given later (table 3). 
Curves for the subject who was observed during two periods (1921 and 
1928) are marked A and B in Figure 2. Both show less curvature than 
all others (c = —.16 and —.20 as compared to —.29 to —.62 for 
others), indicating that such characteristics may persist over long 
periods. 
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FIG. 2—REGRESSIONS ESTIMATED FOR 10 ADULTS INDIVIDUALLY. RANGES AND 
MEANS OF THE REGRESSION COEFFICIENTS ARE INDICATED IN TABLE 1. VER- 
TICAL POSITIONS HAVE BEEN ADJUSTED TO BLOOD UREA = .15 MGM. N PER ML., 
HEIGHT = 170 CMS. DATA FROM VAN SLYKE ET AL (1921, 1928). 


Moller et al report 2, 3 or (for one subject) 5 observations on each of 
nine other subjects. For prisoners of war (15 “normals”, 50 “patients”, 
classified as indicated above) we have two observations of U and V in 
successive hour periods, but only one observation of B which can be 
assumed constant for both hours. No internal regression on b can 
therefore be evaluated for the POW; and for neither group can we fit 
a regression with respectively four and three coefficients for each indi- 
vidual. But we can get an average regression, as in ordinary applica- 
tions of analysis of covariance, from pooled sums of squares and products 
within subjects. Deviations from the average curve will, in these cases 
represent variation between persons. Even when there are more than 
four observations per person it would seem preferable to evaluate the 
average regression for a group in this way rather than by averaging 
individual coefficients; because when observations of independent vari- 
ates are differently distributed there exist variable correlations between 
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the coefficients for each, and because slopes on v in the varying ranges 
contribute information on the quadratic coefficient (the only source of 
an estimate of curvature for an intraclass regression with just two ob- 
servations per class). 

Average regressions thus determined are presented in table 2 for 
five groups of adults, viz. observations made in each of three American 
laboratories and the two groups of POW. Variation is summarised in 
table 3. Mean squares for variation of regressions among individuals 
within each of three groups (subdivision of row 3) are not detectably 
heterogeneous; subdivision of row 4 is discussed above. Row 2 com- 
pared to row 3 indicates heterogeneity of average regressions of the 
five groups, significant at the 5 percent level. However, considering 
the extreme differences of environment, with POW subject to extreme 
protein deficiency, the very small degree of heterogeneity is remarkable; 
furthermore most of it is ascribable to the single individual observed 
by Addis and not to contrast between Americans and POW, (table 2). 
(The latter contrast does not involve regression on b.) Row 3 versus 
row 4 indicates significant heterogeneity of regressions for individuals 


TABLE 2 
No. of | No. of m m' : 
Data Individ- | Observa- at at c k 
uals tions v=Olv=3 
Me er 1 29 —.401 | —.743 | —.570 649 
2. Mclean. ..... 1 17 —.445 | —.650 | —.341 .952 
3. Van Slyke et al: adults} 17 144 — .469 | —.694 | —.375 .847 
4. POW: “normals” . . 15 30 — .555 | —.646 | —.152 _— 
5. POW: “patients” . . 50 100 —.488 | —.640 | —.253 _— 
6. Average all adults 
84 320 —.489 | —.695 | —.344 .869 
7. McIntosh et al: 
. 8 48 —.785 | —.783*| —.017 . 299 
8. Austin ef al: dogs . . 2 26 — .463 | —.826 | —.605 .746 
Standard errors: Rowl ........-. .059 .032 .073 .146 
088 .129 .279 130 
.021 O11 .024 039 


*Linear regression on v, the quadratic term being negligible. 
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TABLE 3 
d.t. M. Sq. 

1. Mean regression (within all adults) . a _ 3 6.44363 
2. Heterogeneity of regressions of 5 groups (table 2) eae 10 .01181 
3. Heterogeneity of regressions among persons within 

groupe... . 96 .005050 
4. Deviations from regressions within persons. asics 127 .002057 


Sub-division of row 3: 


Sub-division of row 4: 


person by 2/8). 25 .002333 
McLean: 1 person (weighted by 1/2) ........2...-. 13 -002234 


*Contains 1 d.f. which properly belongs in row 4, being a deviation from regression for the subject 
with 5 observations, whose regression was not separately evaluated. 


observed in the same laboratory (P < .001), but again the actual mag- 
nitude of differences is small. 

Children observed by McIntosh et al (1928) show similar variation 
(table 5). Their average regression on v appears approximately linear 
with slope corresponding to that for adults at about V = 2.6 ml. per 
minute. This suggests that there should be progressive alteration of m 
with increasing age or size. No such correlation could be demonstrated; 
but, with only 5 to 10 observations per child, individual coefficients 
were erratic. Their regression with blood urea is low: they were not 
given supplementary meat or urea. 

Dogs (Austin et al, 1921) show a similar regression with no coefficient 
significantly different from those for adults. The experiment on dogs 
was run to test correlation of urine urea with blood urea as temporarily 
varied by consumption of urea or meat. For both dogs and men this 
regression is significantly less than unity, the value which has been 
postulated by previous workers. (Table 1 shows that the mean dis- 
crepancy appears still significant when tested against an error based 
on variation between individuals.) 

If it be desired to estimate from an observation of urine urea con- 
centration what it would have been at some other rate of urine excre- 
tion, say V, , most of any required adjustment will be determined by 
the linear coefficient on v. Although regressions exhibited by different 
persons vary significantly, table 1 suggests that the linear coefficient 
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for variation around average v is moderately consistent. Consequently, 
for any observed V not widely divergent from V, , adjustments ap- 
propriate for different persons should be not dissimilar, and may, for 
most purposes, be determined from a well established average curve. 
If there are only one or two observations for a given person, this will 
be the only method of evaluation available. Even with more observa- 
tions, unless unusually numerous, evaluation of an ad hoc regression 
would be unduly laborious and often less accurate to boot (particularly 
if V> should lie outside the range of observations). Regressions on blood 
urea appear somewhat more variable; but their average appears stable 
and normal variation of B is relatively small, therefore similar argu- 
ments may apply. We therefore give in table 2 the average regression 
for all adults. Being derived from observations on both Americans and 
POW’S this regression would appear to have a wide inductive basis. 
Nevertheless it should be noted that neither group is a satisfactory 
sample of a normal population. Therefore although the regression may 
be considered well established and may serve immediate purposes, if 
extensive use is to be made of it it should be re-evaluated for some 
sample properly representative of a population to which it is to be 
applied. 

Using the regression coefficients of row 6 of table 2, urine urea con- 
centrations at a fixed rate of urine excretion (V,) and at blood urea 
concentration (B)) may be estimated from 


log U = log U + .49(log V — log Vo) + .34{(log V)? — (log V,)?} 


— .87(log B — log B,) 
Where U, V, B are observed values. If the adjustment for B is included 
it should be noted that this will estimate urine urea only for B changed 
to B, by a temporary alteration in diet; the regression does not hold for 
long term alterations. If it is required to adjust a mean of several 
observations, means of logarithms should be used and the value in- 
serted for (log V)? should be the mean of (log V)’ and not the (mean 
of log V)’. 
If comparison is to be made between persons of different sizes a 
further term, —1.5 (log H — log H)), may be added, where H = height 
(see below). 


EXTERNAL REGRESSIONS 


By a sub-division of analyses of variance, similar to that given by 
Kendall (1948, Sec. 24.30), we may examine whether variation between 
men indicates any different relationship to the independent variates as 
compared to that within. This analysis is given in table 4. For van 
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Slyke’s adults we consider heterogeneity of regression on the three 
variates (v, v’ and b); POW’s, for the reason given above, show internal 
regression only on v and v’. Besides these table 4 shows the part of 
sums of squares between men attributable to regression on an extra 
variate for which there are only single observations per man; for van 
Slyke’s adults on height, for POW’s on blood urea. (Heights of POW’s 
were not recorded.) 


TABLE 4 
v. Slyke et al: Prisoners of war 
adults normals patients 


df. | M.Sq.| df. | M.Sq.| df. | M.Sq. 


Heterogeneity of regres- 
sions “within” and 


“between” ..... 3 .02007 2 .00095 2 .04295 
Subsidiary regression 

(see text) ..... 1 .02408 1 .04413 1 .04310 
Deviations from regres- 

sion between men . . 12 .02221 11 .01817 46 .04280 


Total between men (“re- 
duced” for adjustment 
by internal regression) 16 .02193 14 .01756 49 .04281 

Deviations from average 
internal regression . 124 .00243 13 .00263 48 .00672 


No different relationship between as compared to within men can be 
detected. The enhanced mean square for adjusted means apparently 
represents individual variation occurring at random relative to the in- 
dependent variates investigated. 

Between POW’s no correlation with blood urea concentration can 
be detected. (The condition that the mean square for this regression in 
the group of normals is greater than the remainder mean square—albeit 
not significantly so—is due entirely to a single observation. For both 
groups together the correlation would be approximately zero. See 
figure elsewhere—Smith, 1949.) For van Slyke’s adults this regression 
is similar both between men (k = .837 + .109) and within (.847 + .039). 
Since reasons have been given in the earlier paper for believing that the 
normal relationship between individuals on unequal stable diets may be 
as in the POW’s, the relationship shown by van Slyke’s data calls for 
comment. Major differences between the blood ureas of his subjects 
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were due to some having consumed urea, and consequently reflect the 
same temporary fluctuations which determine the internal regressions. 
Among means of thirteen persons who were not given urea the regres- 
sion, after adjustment for V, is k = .38 + .22. It seems likely that 
persons on normal diet may show some, but only low, correlation between 
blood and urine urea concentrations. If it be still considered worth 
while to compare urine urea concentrations this point requires further 
investigation on a larger sample which is properly representative of a 
definable population. 


CORRELATION OF UREA EXCRETION WITH BODY SIZE 


For most adults the effect of body size on urine urea concentration 
at fixed V is usually negligible relative to random variation. But for 
persons divergent from normal (taken as 170 cm. in height) adjustment 
for size has been recommended. Méller et al (1928) presented reasons 
for considering that the best characteristic for this purpose is surface 
area, which they computed from height and weight by the formula of 
Du Bois. Later McIntosh et al (1928) argued that weights of patients 
suffering from disease are unreliable as an indication of their normal 
sizes, and they prefer to base size adjustments on height alone. Their 
procedure is, however, circuitous and confusing. 

From height they estimate what they term the body weight “ideal’’ 
for height, by which, apparently, they mean average body weight for a 
population of given height. The estimate of average body weight thus 
obtained is then recombined with the height measurement to estimate 
surface area. Actually the only variate is height, and the net effect is 
the same as if adjustment were based directly on it alone. To introduce 
a supposititious surface area is merely to mask what is in fact the basis 
for the corrections, and to complicate computations. 

To investigate whether weight, either in its own right or as an indi- 
cator of surface area, may add to the accuracy of adjustment, more 
information may be gained by considering both variates (height and 
weight) explicitly, rather than arbitrarily combined in a formula under 
another name. The experiments on children (McIntosh et al, 1928) 
were planned specifically to consider effect of body size. After adjusting 
means per child by their internal regression on v and b, regression on 
size characters gives analysis of variance as in table 5. This shows that 
once adjustment has been made for height, nothing further is gained by 
considering weight. (Correlation with log weight alone is nearly, but 
not quite, as great as with log height. It would account for .708 in the 
sum of squares, with regression coefficient 0.97 + .24.) 

Variation in size of the 17 adults observed by van Slyke ef al is too 
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TABLE 5 
d.f. S. Sq. M. Sq. 
Between children (adjusted for v and b) . . . 7 .9691 
Regression on log height. . . ...... 1 . 7881 
Deviations from regression on height . . . 6 .1810 |° .03018 
Partial regression on log weight ..... 1 .0107 
Within children 
Heterogeneity of internal regressions .. . 21 .00712 
Deviations from individual regressions . . . 3 .00274 


small to show significant correlation (subsidiary regression, table 4). 
However, the regression is not dissimilar to that observed for children, 
and here also correlation with height was slightly. greater than with 


weight. 
Regression coefficients of u on log height are: 
for children 1.66 + .33 
for adults .96 + .92 
weighted mean 1.59 


Since surface area as computed by van Slyke et al appears to be pro- 
portional to (height)'**° for children, and to (height)’*** for adults the 
hypothesis on which their adjustment is based appears reasonable. 


POSTSCRIPT 


A referee, on reading an earlier version of this paper, suggested 
alternative formulation which appears of interest and worth comment. 
He suggests that the important quantities to be considered are: 


total urea excreted per minute = UV 

concentration potential = log U — log B, 

volume rate = V, 

and normal liquid handling capacity of the kidney tubules: 


and he asks, very pertinently, why was U chosen as the significant de- 
pendent variate? 

Relative to the van Slyke test the total amount, UV, is more im- 
portant than the concentration U. But V being still the independent 
variate, assumed for purposes of curve fitting free of error, to use either 
u or (u + v) as dependent variate leads to identical regressions on v 
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(except for a difference of unity in the linear coefficient). Consequently 
it is simplest to deal with the variates observed. According to con- 
venience in any particular context the regression may still be expressed 
in either form. 

Secondly he suggests that if ratio of volume excreted to normal 
capacity may be important, and if we may assume normal capacity « 
weight « height*, we should use as independent variate V/H*. Else- 
where (1949) I have criticised the a priori assumption that the ratio 
U/B is significant, and (1936) pointed out that to use arbitrary ratios 
may usually be inefficient. The first need it to establish the relative 
effects of both variates. If confidence intervals of regression coefficients 
include values conforming to some combination of the variates which 
has theoretical backing, then the combination may be adopted if it 
gives a simplification. For investigation one could use (V/H*) and 
either V or H as the two variates. But it seems simpler initially to use 
the primary variates actually observed, rather than to focus attention 
on a particular one‘of many possible hypotheses. We may then consider 
what would be the effect on the straight-forward regressions if some 
compound of the variates were in order. Linear combinations present 
no problem since corresponding terms of the regression can be freely 
re-arranged as in the preceding paragraph. To examine the effect of 
introducing a quadratic term such as c(v — ah)’—where h = log H, a 
is some constant and we suppose c ~ 0—, we would note that, if this 
may be a better fitting regression than (6), then m (or m’) for each 
individual should show regression on h with slope —2ac. Secondly the 
partial linear regression on h when evaluated in the usual way should 
be —am’, where m’ = du/dv near the mean of v, for our data about —.7. 
Therefore, if (V/H*) is to be used as independent variate, both adults 
and children indicate that a should be about 2. However variation of 
H among adults is too small to show any effect of incorporating H in 
the quadratic term, and for children c is approximately zero so that the 
whole term drops out anyway. Therefore no evidence as to the merits 
of this hypothesis can be deduced from present data. 


SUMMARY 


Urea clearance tests based on the van Slyke formulae were found to 
be very misleading when used in a prisoner of war camp where diet was 
deficient in protein. Subsequent comparison of the POW observations 
with others recorded in the literature revealed that the test appears to 
have been based on (statistically) false premises. Before this was 
realized considerable statistical analyses had been done, both of the 
POW data and of those on which the test was founded, with a view to 
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simplify the test. Although that purpose is now irrelevant an account 
of the work is given here in the hope that it may aid further research. 

The relation between urea concentration in urine and volume of 
urine excreted in a given time varies significantly between one person 
and another. All, however, show similar curves which can be approxi- 
mated by a quadratic in logarithms of the variates. 

The regression of log urine urea on log blood urea concentration is 
about .85 for temporary fluctuations. This is significantly less than 
unity as postulated by previous workers. Between persons on perma- 
nently different diets these variates may be uncorrelated. 

Urine urea concentrations of children at fixed urine excretion rate 
vary approximately as (height)*”’, or in proportion to surface area as 
postulated by previous workers. 
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INTRODUCTION 


ALCULATION of the fiducial limits of error of a bioassay is intricate. 

A technical assistant, computing by rote, is unlikely to be able to 

decide whether an abnormal result is due to miscalculation or is in- 

herent in the data. Nomography and the use of error curves can reduce 

his labors, promote his confidence, and advance his insight into the 
factors contributing to the assay error. 

Several graphical aids for bioassays and errors are already available, 
e.g. those of de Beer [1], Knudsen & Randall [2], Bliss [3], Sherwood 
[4], Goyan & Dufrenoy [5], Koch [6], and Harte [7]. Lees’s specific 
calculating machine [8] should also be mentioned in this connexion. 
The most recent aid, and perhaps the one most nearly suited to the 
requirements of the writer’s laboratory, is that of Healy [9], whose 
nomogram gives errors for assays of given size and structure. The 
present requirements, however, were error charts covering 4- and 6- 
point assays of given design but greatly varying size. It was not re- 
quired to calculate graphically the mean estimate itself, the writer 
being unable to appreciate the advantages of so replacing the simple 
arithmetic operations involved. 

The devisal of appropriate error charts, the main subject of this 
paper, has thrown into relief some aspects of assay design and error 
factors that may be worthy of concomitant record. 


ERROR THEORY 


To explore the graphical possibilities of the error distribution, let 


“us derive a comparatively simple general expression for the fiducial 


*Presented at a joint meeting of the Biometric Society (British Region) and the Biological 
Methods Group of the Society of Public Analysts, London, March 14, 1950. 
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RESPONSE 


LOG DOSE [RELATIVE] 


FIGURE 1. 
Showing the geometric relations between J, D and B for a 4- or 6-point (or any 2k-point) assay. Either 
set of linear points may represent the standard preparation (S), the other then representing the test 
material (U). 


limits of a continuous response assay with equal numbers of responses 
in each of 2k dosage groups, where k is any integer other than unity. 
Half of the groups shall be on the standard preparation and half on 
the test material, and there must be a uniform potency ratio, antilog I, 
between the doses for adjacent groups of either preparation. Let n be 
the total number of responses and s the pooled estimate of the homo- 
geneous standard deviation of responses within groups. Assume a linear 
relation between response and log dose. Let B be the mean response 
increment between adjacent dosage groups of either preparation (i.e. 
the coefficient of regression, or slope, b, is B/Z), and let D be the corre- 
sponding mean increment of all ‘test’? responses over all “‘standard”’ 
responses. The fraction D/B, which we shall make = A, determines 
the best estimate R of the unknown potency of the test material in 
terms of the standard (dose for dose); for R = antilog M = antilog JA. 
Figure 1 illustrates these-relations. 

Now A can be regarded as a sample of the population value a be- 
longing to an infinitude of assays identical in size, design etc., with the 
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actual one. We can restate this as the long-run tendency for aB to 
be equal to D, and this consideration provides the clue, first followed 
up by Fieller [10], to the required fiducial limits of A and, consequently, 
of R. The difference D — aB is normally distributed about a mean of 
zero and therefore, when measured in units of its estimated standard 
deviation, will follow the ¢ distribution with the number of degrees of 
freedom associated with s. The estimate of the standard deviation of 
D — aBis 


~/V(D) — 2aCov(DB) + V(B), 
so that we may write 
t? = (D — aB)’/[V(D) — 2aCov(DB) + a’ V(B)], 
which is a quadratic in a: 
a’[B’ — ?V(B)] — 2a[BD — Cov(BD)] + — #V(D) = 0. 
This equation can be simplified if we introduce Fieller’s statistic, 
C = B’/(B’ — #V(B)), 


and make use of the facts that, in assays of the type we are here con- 
cerned with, 


V(D) = 4s°/n, V(B) = 12s’/n(k’ — 1), 
and the covariance of B and D nil. It follows that 
— ?V(B) = B’/C, 
and that 
?V(D) = — — 1)/3C. 
Substituting these in the above equation, multiplying throughout by 
C/B*, and replacing D/B by A, we finally arrive at 
— 2aCA + CA? — (C — 1)(k — 1)/3 = 0, 
whose roots are 
CA + — 1][CA® + — 1)/3]. 


The roots give the upper and lower fiducial limits, i.e. the range 
whose probability of not enclosing a is P, a value determined in the 
choice of ¢ and usually taken at 0.05. To obtain the percentage fiducial 
limits we must use 


ACC 1) 2 VC — + — 1/3), (1) 
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which gives two values that, multiplied by J and added to 2, will be 
the logs of the required percentages. 

Several points of interest attached to this expression are worth 
special mention: 
(i) When A is zero, i.e. when the mean responses to the standard 
preparation and the test material are equal, it reduces to 


and illustrates the biometrical advantage of shrewd guessing at the 
potency of the test material, and consequently of so choosing the 
dosages to bring A as close to zero as possible. It is also apparent that 
only when A is zero will the fiducial limits of the logarithm of the 
assay be symmetrical about the mean; as A moves away from zero a 
progressive bias sets in. Note that the expression takes the very simple 
form + +~/(C — 1) for a 4-point assay with equal responses. 
Theoretically A has no upper or lower limit, but in practice it is 
seldom that A can exceed +2 without straining the validity of the 
assay, usually by a breakdown of the parallelism of the “standard” and 
“test”’ response curves. 
(ii) Because in theory the variance of the slope of the mean response 
curve can have any value including 0, as (’V(B)/B’ passes from zero 
in a positive direction, so C passes from unity to + and then from 
— to zero, following the type curve y = 1/(1 — x). As C approaches 
unity the value of the expression dwindles towards vanishing point 
where the precision of the assay apotheosizes at 100%. A value of 
infinity, or even a negative value, for C is, by contrast, realizable. At 
infinity the position is interesting: it is reached when B’ = (’V(B), ie. 
when the slope of the mean response curve is exactly poised between 
significant and non-significant differences from zero at the selected 
probability level. It is then obvious that the assay is almost worthless, 
and this is reflected in the shift of the bias (the mid-point of the error 
range) to infinity, and the expansion of the range itself to nearly + ©; 
in fact one of the two values of our expression becomes infinity and the 
other a comparatively small negative whose magnitude depends on A 
and k. This latter value can readily be found if we equate expression 
(1) to Z and rearrange the equation so that C is collected. We have 


1-—1/C = L?/[((L + A)? + — 1/3); 


whence, if C = ~, we obtain 


L = —[A? + (k? — 1)/3]/2A 
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for the required value of Z. Remembering that the other value is 
infinity (positive or negative according to the sign of A), and that we 
have to take antilogs in the transformation of Z to limits of error, we 
can see that when A is positive (i.e. when the test material is more 
potent than hypothesized) the upper fiducial limit will be infinity but 
the lower limit will have a finite value. On the other hand, when A 
is negative (i.e. when the estimated potency is low) the lower limit will 
be nil but the upper limit finite. Only in the special case of A = 0 
will the range extend from 0 to . Reflection on this critical state of 
slope significance will show that the emergence of one finite limit is 
eminently reasonable. If, for example, a material expected to contain 
100 units per gram assays at 200 units per gram, with a slope of enig- 
matic significance at P = 0.05, we cannot expect to extract any in- 
formation on how large the true unitage may be, yet we can be fairly 
confident that the material has some potency, and we shail not be sur- 
prised to find that, in fact, a unitage is calculable below which there is 
only a 2.5% probability of encountering the true value. 

Occasionally C may be negative; this means that V(B) > B’ and 
that there are no real fiducial limits. (The equations for Z are soluble 
with negative C’s under certain conditions, but the roots so obtained 
are imaginary.) The assay is consequently entirely valueless at the 
working probability level. 

(iii) The two commonest values for k are 2 and 3, so we have: . 


for the 4-point design: A(C — 1) + VW/(C — 1)(CA’ + 1); 


for the 6-point design: A(C — 1) + W/(C — 1)(CA’ + 8/3). 


The comparatively simple relation of these two expressions masks 
‘complex and contingent relations between the errors of the two designs 
for the same number of subjects (animals). One simple deduction can 
however be made: for given values of n, s and b a 6-point design with 
interval J tends to biometrical identity (interactions or validity tests 
aside) with a 4-point design with interval I +/8/3 for all values of A; 
the tendency is a function of ¢, whose degrees of freedom will differ 
slightly for the two designs, the difference becoming negligible as n in- 
creases. In general the 4-point assay is more or less efficient than the 
comparable 6-point assay if the interval is > or < I +~/8/3. This 
aspect of design is enlarged on in the Appendix. 
(iv) Finally, if, as above, we let L denote a root of the expression, it 
is manifest that error curves applicable to given values of LZ can be 
drawn with C and A, or with functions of C and A, as co-ordinates. 
This is the basis of the charts discussed hereunder. 
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BACKGROUND OF CHARTS 


Of the two co-ordinates for the projected L curves, A is a simple 
key datum that is normally calculated en route to the evaluation of 
the assay itself. It may, in fact, be preferable to graduate the A scale 
as antilog JA, which is R, the best estimate of the potency of the test 
preparation as a multiple of the expected potency (the null hypothesis), 
and this practice has been adopted here. On the other hand C, which 
is B’/[B’ — 12¢’s’/n(k’? — 1)], is complex, a function of five variables. 
However, one of these, ¢, is wholly determined by n and the structure 
of the assay, and another k, can be fixed; this leaves three variables: 
B, s and n. Now if we take the trouble to calculate an appropriate 
function of any one pair of these, we shall be left with two values to 
derive C from, and this lends itself to simple nomography. The best 
pair for choice is undoubtedly B and s, because they are components 
of the important statistic 4, whose evaluation is commendable per se. 
Lambda is the residual standard deviation of the response measured in 
slope units i.e. 


= s/b = 8 I/B, 


and being a ratio, independent of the kind of response, of the size of 
the metameter, and of n (the size of the assay), \ is a coefficient of the 
inherent sensitivity of the particular set of animals used in an assay. 
To evaluate \ as a matter of routine, to keep track of its fluctuations 
over time, and to compare it with evaluations in other laboratories, is 
a useful discipline. To know the relative mean \’s of two types of assay 
conducted with the same kind of animal is to know the relative cost, 
as far as animal usage is concerned, of the two assays at any degree of 
precision. And the two components of \ are necessarily always avail- 
able. Bliss and Cattell [11] have compiled a useful list of \’s for various 
kinds and specimens of bioassay. 
Introducing \, we may re-fashion the equation for C as, 


C = 1/fi — 12°? — 1), 


showing the basis of its derivation from two variables \ and (’/n. In 
sum, then, for particular values of J and k, we can construct a nomogram 
to estimate C from \ and n, and use the C scale as an axis of the error 
curves. Thus for an assay of given design the percentage fiducial 
limits of error at a given probability level can be read off from the three 
data, \, n and A. It may be well to emphasize that the chart as an 
entity, nomogram plus graphs, must have J and k determinate—we 
cannot, for instance, use one set of error curves in conjunction with 
two nomograms based on different J’s and/or different k’s. 
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CHART OF PERCENTAGE LIMITS OF ERROR (P+0-05) OF 4-POINT BIO-ASSAYS 
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The two specimen charts herewith reproduced cover the 4-point 
assay with 1:2 “steps” (k = 2; 7 = 0.30103) and the 6-point assay with 
3:5 “steps” (k = 3; J = 0.22185). The centre horizontal scales of F(C) 
are graduated arbitrarily as in practice they act merely as gate-ways to 
the error curves. Choice of a suitable function of C turns on several 
considerations including the desirability of a linear distribution and the 
accommodation of suitable ranges and spacings of )’s, n’s and the error 
curves. The best choice is ~/(C — 1)/C, which is the square root of 
Finney’s statistic g [10]. 

The nomogram for the 4-point assays has a double n scale: one for 
“no blocks’, i.e. ¢ is based on n — 4 degrees of freedom; the other for 
blocks of four (e.g. isogenic tetrads of subjects), t being based on 3n/4 — 
3 degrees of freedom. As blocks of 6 are rare in bioassay work, only 
the blockless design has been catered for on the 6-point chart. 

Symmetry of design, i.e. equal numbers of responses at each point, 
is a theoretical condition of the use of either chart; nevertheless the 
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CHART OF PERCENTAGE LIMITS OF ERROR (P=0:05) OF 6-POINT BIO-ASSAYS 
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slight asymmetry brought about by a few lost responses will not seriously 
distort the error from that that would have been obtained by the same 
losses evenly distributed among the groups. Indeed, to judge from some 
examples given later in this paper, asymmetry must be considerable 
before the charts cannot safely be used to find the fiducial limits. 


CONSTRUCTION OF CHARTS 


The nomogram consists of two parallel uniform scales connected by 
a third intercept scale to form an N-shaped triad. If the parallel scales 
are zero’d at the junctions, the diagonal scale can be calibrated to yield 
ratios of aligned pairs of readings. For the present purpose it is meet 
to distribute \ on one parallel and +/g on the other, and from them to 
calibrate the diagonal in n (necessarily non-uniform). Then a straight 
edge laid across the diagram and connecting any \ with any n will cut 
the +/g scale at the correct value. 

The size of the scaling must first be settled. Experience of the 
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order of magnitude of \ likely to be met with, together with considera- 
tions of the lay-out of the chart as a whole, indicate that a \ scale from 
0 to 0.5, coupled with a 4/9 scale of 1 to 0, is suitable. It should be 
_ noted that g = +/g = 0 corresponds to C = 1 and absolute precision, 
that g = +/g = 1 corresponds to C == and only one real fiducial 
limit, and that g > 1 results in a negative value for C and absolute 
imprecision. A perpendicular drawn between these two scales will cut 
and +/g where the latter is numerically equal to 1 — 24. Now as 


= — 1)/12, 


the point on the diagonal cut by the perpendicular must correspond to 
a unique value of n for a given J and k. This provides the key to the 
calibration of the diagonal. Replacing V9 by 1 — 2d in the above 
expression, and re-arranging, we find 


— 1)/12/[t + — 1)/12], 


from which any required n can be located on the diagonal as the inter- 
cept of the perpendicular from the corresponding ». To locate, for 
instance, n’s in the “no block” scale for J = log 2 and k = 2 at the 
5% probability level we substitute in 


= 0.30103 ~/n/(2t + 0.60206 ~/n), 


the ¢-table being entered at n — 4. Thus for n = 20 the value of X is 
0.1942; so we locate this n as the intercept from the perpendicular at 
this \. In practice it is necessary to locate only a few n’s by calculation 
and construction, the others being found by graphical interpolation. 

The +/g scale must now be considered as one of the ordinates for 
the error curves to be drawn alongside. As the 4/9 scale thus acts 
merely as a stepping stone from the nomogram to the curves, in the 
finished diagram it need be neither labelled nor figured. The other 
ordinate shown on the finished charts is R, which is “found potency of 
test material as a multiple of the hypothesized potency”, i.e. antilog 
TA, and the range chosen is 0.5 to 2, which will cover the great majority 
of assays. (In practice assays yielding values outside this range are 
generally repeated with an amended hypothesized potency.) The 
working scale however is the more amenable A, distributed linearly’ 
between +1 if I = log 2, and between +1.357 if I = log 5/3. 

We require to plot A against ~/g for given values of 


A(C — 1) 4+ V(E — + — 19/3] 


— or, more particularly, for given percentage limits of error at what- 
ever probability level has been built into that expression and which 
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here happens to be P = 0.05. Let L be a root of the expression. Re- 
calling that C = 1/(1 — g) we can easily derive that 


Vg = +L/V(L+ AP + 1/3. 


Let us consider the curve of all combinations of A and +/ g that 
yield an upper limit of 160%, i.e. 1.6 times the mean assay estimate, 
for 4-point designs (k = 2) and J = log 2. The logarithm of 1.6, which 
is 0.2041, is JL. Therefore L = 0.6780. We can now calculate VJ g 
for any value of A. For example: 


A v9 


1.0 0.3471 
0.5 0.4388 
0 0.5612 
—0.5 0.6675 
~19 0.6454 


The whole curve can be drawn from these five points. It is use- 
ful however to include a sixth, viz, the maximum, which always 
occurs at L = —A, and which for this curve is at A = —0.678 and 
/g = 0.678. Incidentally, the maxima of all the curves fall on a 
straight line drawn from the right centre (A = 0; /g = 0) to the 
points A = + ~V/(k* — 1)/3 and ~/g = 1, which are the corners of 
the frame of the curves for assays with dosage intervals of 1:2, and 
k = 2, but not otherwise. 

Alignment of the maxima is one of several devices that can facilitate 
the construction, and make the task less formidable than it may super- 
ficially appear. For example, the points for particular curves need not 
be calculated in sets as outlined above; it is best to begin with a tabula- 
tion of L and then to calculate ~/g for all points on the line A = 1, 
then for A = 0.5 and so on. Nor need “all points” be taken literally: 
interpolation can be used to lighten the arithmetic. (It is nevertheless 
advisable to construct a few “tracer” curves on a good number of 
points; this helps accuracy of flow of the other curves.) Another helpful 
circumstance lies in the fact that the pattern of the curves of upper 
limits is a mirror image of that of the curves of the lower limits, and 
many individual pairs of curves are mirrored; for example 50% and 
200%, 82% and 122%. 


VALIDITY OF ASSAY DATA 


A pre-requisite for the use of the charts, and indeed for any method 
of estimating the mean potency and its precision, is that the disposition 
of the responses at the 2k points shall be consistent with the basic 
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order of magnitude of \ likely to be met with, together with considera- 
tions of the lay-out of the chart as a whole, indicate that a \ scale from 
0 to 0.5, coupled with a +/q scale of 1 to 0, is suitable. It should be 
_ noted that g = 4/9 = 0 corresponds to C = 1 and absolute precision, 
that g = +/g = 1 corresponds to C = and only one real fiducial 
limit, and that g > 1 results in a negative value for C and absolute 
imprecision. A perpendicular drawn between these two scales will cut 
hand +/g where the latter is numerically equal to 1 — 24. Now as 


V9/d = t/IV — 1)/12, 


the point on the diagonal cut by the perpendicular must correspond to 
a unique value of n for a given J and k. This provides the key to the 
calibration of the diagonal. Replacing V9 by 1 — 2A in the above 
expression, and re-arranging, we find 


= — 1)/12/[t + — 1)/12], 


from which any required n can be located on the diagonal as the inter- 
cept of the perpendicular from the corresponding ». To locate, for 
instance, n’s in the “no block” scale for J = log 2 and k = 2 at the 
5% probability level we substitute in 


= 0.30103 ~/n/(2t + 0.60206 /n), 


the ¢-table being entered at n — 4. Thus for n = 20 the value of d is 
0.1942; so we locate this n as the intercept from the perpendicular at 
this \. In practice it is necessary to locate only a few n’s by calculation 
and construction, the others being found by graphical interpolation. 

The +/g scale must now be considered as one of the ordinates for 
the error curves to be drawn alongside. As the V9 scale thus acts 
merely as a stepping stone from the nomogram to the curves, in the 
finished diagram it need be neither labelled nor figured. The other 
ordinate shown on the finished charts is R, which is “found potency of 
test material as a multiple of the hypothesized potency”, i.e. antilog 
IA, and the range chosen is 0.5 to 2, which will cover the great majority 
of assays. (In practice assays yielding values outside this range are 
generally repeated with an amended hypothesized potency.) The 
working scale however is the more amenable A, distributed linearly: 
between +1 if J = log 2, and between +1.357 if I = log 5/3. 

We require to plot A against ~/g for given values of 


1) 4 V(C — 1[CA’ + — 1)/3] 


— or, more particularly, for given percentage limits of error at what- 
ever probability level has been built into that expression and which 
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here happens to be P = 0.05. Let LZ be a root of the expression. Re- 
calling that C = 1/(1 — g) we can easily derive that 


Vg = +L/VL+ A) + 17. 


Let us consider the curve of all combinations of A and +/g that 
yield an upper limit of 160%, i.e. 1.6 times the mean assay estimate, 
for 4-point designs (k = 2) and J = log 2. The logarithm of 1.6, which 
is 0.2041, is IL. Therefore L = 0.6780. We can now calculate +/g 
for any value of A. For example: 


A V9 


1.0 0.3471 
0.5 0.4388 
0 0.5612 
—0.5 0.6675 
0.6454 


The whole curve can be drawn from these five points. It is use- 
ful however to include a sixth, viz, the maximum, which always 
occurs at L = —A, and which for this curve is at A = —0.678 and 
~/g = 0.678. Incidentally, the maxima of all the curves fall on a 
straight line drawn from the right centre (A = 0; /g = 0) to the 
points A = + +~/(k’ — 1)/3 and ~/g = 1, which are the corners of 
the frame of the curves for assays with dosage intervals of 1:2, and 
k = 2, but not otherwise. 

Alignment of the maxima is one of several devices that can facilitate 
the construction, and make the task less formidable than it may super- 
ficially appear. For example, the points for particular curves need not 
be calculated in sets as outlined above; it is best to begin with a tabula- 
tion of L and then to calculate ~/g for all points on the line A = 1, 
then for A = 0.5 and so on. Nor need “all points” be taken literally: 
interpolation can be used to lighten the arithmetic. (It is nevertheless 
advisable to construct a few “tracer” curves on a good number of 
points; this helps accuracy of flow of the other curves.) Another helpful 
circumstance lies in the fact that the pattern of the curves of upper 
limits is a mirror image of that of the curves of the lower limits, and 
many individual pairs of curves are mirrored; for example 50% and 
200%, 82% and 122%. 


VALIDITY OF ASSAY DATA 


A pre-requisite for the use of the charts, and indeed for any method 
of estimating the mean potency and its precision, is that the disposition 
of the responses at the 2k points shall be consistent with the basic 


we. 
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assumptions of the assay [13]. Validity tests must therefore be applied. 
(Strictly, they test for non-validity; validity is unprovable, but is 
assumed in the absence of evidence to the contrary.) For the 4-point 
assay there is only one test available: that for dissimilarity of the slopes 
of the response curves of the standard preparation (S) and test material 
(U). If we call the four response means S, , S, , U, and U, , the fraction 


Si S. + U; U2 
2s/+/n 


must, if the assay is unequivocally acceptable, be not more than the 
P = 0.05 value of ¢ in Fisher’s table, suitably entered. If it is more, 
the assay may or may not be acceptable; the choice turns on whether 
or not there is extraneous information on the characteristics of the re- 
sponse curves; and to expand this statement we must first consider 
the validity tests for the 6-point design. 

Three distinct tests can be applied to a 6-point assay. Firstly, as 
there are three points per curve, the assumption of linearity of the 
combined S and U curves can be tested: non-linearity should entail a 
significant difference between the upper and lower “‘steps”. Secondly, 
the assumed parallelism of the two curves can be tested: non-parallelism 
should entail a significant difference between the distances apart at the 
ends, i.e. between U, — S, and U; — S;. Thirdly, the shapes, assumed 
to be identical, of the two curves can be compared: contrasting curvature 
should show up as a significant difference between, on the one hand, the 
distance apart at the middle and, on the other hand, the mean of the 
distances apart at the ends, i.e. between U, — S, and (U, + U;)/2 — 
(S, + S;)/2. So the three fractions that we can apply the ¢ test to 
(with n — 6 degrees of freedom) are: 


numerator denominator 


for linearity: 

(U, — 2U2 + Us) + (Si — 282+ 8-V72/n; 
for parallelism: 

(U, — Us) — — s/24/n; 
for common curvature: 

(U, — 2U, + U3) — (8S, — 282+ /n. 


Illustration of possible dispositions of the 6 points is best done with 
(idealized) curves on an absolute, rather than on the customary relative, 
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RESPONSE 


LOG DOSE [ABSOLUTE] 


FIGURE 2. 


Showing the ideal disposition of the two response curves S and U, and three possible departures from 
the ideal. 


I coincident and linear, 
II coincident and quadratic, 
III non-coincident and linear, 
IV non-coincident with dissimilar curvature, 
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log-dose scale, as in Figure 2; and in consequence “parallelism”? must 
be transposed to “coincidence”. In perfection S and U have a common 
rectilinear response/log-dose curve (1); all validity tests are then passed. 
But S and U may share a common quadratic curve (II); it is obvious 
that there is no fundamental objection to this situation; a finding of 
non-linearity is therefore not a true test of validity, but it may mean 
that a transformation of the metameter will improve the assay. The 
third possibility is that S and U have distinctive rectilinear curves (III); 
then “non-parallelism”’ emerges from the appropriate validity test; this 
means that the standard preparation and the test material are acting 
in biochemically different ways, and that the basis of comparison is 
unacceptable. Finally, one relation may be linear and the other curvi- 
linear (IV), or both may be curvilinear in different ways; this results in 
failure of the “common curvature” test, and it disturbs the basis of 
comparison almost as strongly as non-parallelism. There are, of course, 
combinatorial possibilities, e.g., non-parallelism in association with con- 
trasting curvature (which association, incidentally, is depicted in set 
IV), and they will be reflected in corresponding linkages in the validity 
tests. 

Now the relative positions of S, , U, , S; and U; in sets II, III 
and IV of Figure 2 have deliberately been made identical. Therefore, 
if we ignore S, and U, in these sets we are left with three 4-point assays 
with underlying characteristic differences that would nevertheless be in- 
distinguishable with the data available. Herein lies the ambiguity 
occasioned by failure of the “‘validity” test of a 4-point assay. Failure 
may be accounted for by the existence of a quadratic, instead of a linear, 
response curve, and this is not an invalidating circumstance, for not 
only is it biometrically unexceptionable, but it so happens that disregard 


‘of the quadratic term does not falsify the estimate—the linear equation 


applied to the 4 or 6 points yields the same result as the “true” quadratic 
equation. 

If we can be confident from the inferences of other experiments that 
the S and U curves are coincident for this type of assay, the 4-point 
design may be preferred on the score of simplicity or, in certain cir- 
cumstances, of efficiency (see Appendix). 


OPERATION OF CHARTS 


Given the response data the assayist with an appropriate chart avail- 
able has six steps to take to the full assay: 

(i) Calculation of mean group responses 

(ii) Estimation of standard deviation 

(iii) Testing for assay validity 


| 
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(iv) Calculation of 
(v) Calculation of “found/expected potency of test material” and 
hence of the best estimate of absolute potency in whatever 
metric is customary 
(vi) Entering the error curves through items (iv), (v) and the total 
number of responses, to read off the percentage fiducial limits. 
The calculation of the mean group responses needs no elaboration. 
Let us call the means S, , S, ete and U, , U, ete for, respectively, the 
standard preparation and the test material, at the separate log dosage 
intervals 7. Estimation of the standard deviation of responses within 
groups is also straightforward, although if there are blocks, e.g. litters, 
across the groups we must eliminate the inter-block variance by 
means of a simple analysis of variance. Let the estimate be s. The 
associated number of degrees of freedom will normally be n — 4 for 
4-point assays and n — 6 for 6-point assays, where n is the total number 
of responses, but appropriately fewer if blocks have been taken into 
account, e.g. 3n/4 — 3 for a 4-point assay. 
Having decided the validity question (see previous section), we 
next calculate \: for the 4-point design with dosage steps of 1:2, we have 


0.6021s 
|-U,+ U.— 8S, + 
and for the 6-point design with dosage steps of 3:5 we have 
0.88748 : 
| -—U, + U; — 8, + 8S; | 


The ratio R, labelled ‘“found/expected potency of test material” on 
the charts, is 


U,+ U,— 8, — 


antilog 0.30103 


for the 4-point “1:2” design and 


U,+ U,+ U; — 8, — 8 — 8S; 
—U, + Us — 8, + Ss 


for the 6-point “3:5” design. Multiplication of the expected (null 
hypothesis) potency by the ratio give the best estimate of potency. 
Finally we enter the chart at the appropriate A, cross to the centre 
scale with a straight edge aligning \ and n, and from there move verti- 
cally down into the error curves until the correct horizontal is inter- 
cepted. The upper and lower fiducial limits are read off at that point. 
In proceeding along the above lines we ignore, by implication, any 


antilog 0.2958 


a 
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TABLE 1 
RESPONSES OF A CHICK ASSAY OF VITAMIN D 
Si U; S: Uz S; U; 
115 124 79 121 80 100 
141 163 78 163 88 89 
129 134 104 160 80 85 
130 160 93 123 90 73 
137 112 69 136 104* 92 
112 151 86 139 87 79 
140 154* 109 158 72 85 
158 134 113 109 93 111 
148 149 130 138 95* 101 
147 165 75 122 87 105 
110 144* 107 116 95 112* 
122 167 98 67 75* 89 
142 141 85 102 72 120 
132 149 122 102 86 82* 
142 161 140 99 72 89 
91 150 122 123 66 92 
120 150 116 47 66 124 
148* 144 91 125 93 80* 
129 160 94 108 99 72 
154 150* 93 129 110 84 
167 151 100 121 81 87 
121* 155 110* 93 63 107 
140 153 106 121 103 90 
142 139 103* 117 102* 100 
124 159 119* 90 82 68 
144* 141 102 153 101* 73 
127* 144 95 149* 84 110 
165 130 87 144* 74* 98 
Tora.s 3,777 4,134 2,826 3,375 2,400 2,597 


N.B. For the meaning of the * and of the trisection of the 168 data, see text. 


distortion introduced when n is not evenly sub-divided among groups, 
e.g. when some responses of a symmetrically planned assay are not 
available at the end of the experimental period. Asymmetry will effect 
D/B (the contributory mean response difference should be appropriately 
weighted), and therefore the mean assay estimate, and will effect V(D) 
and V(B) to displace the fiducial limits. In practice however the effects 
are generally too fine to be of importance—see next section. 
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SOME EXAMPLES 


The data from a sizable 6-point chick assay of the vitamin-D 
potency of an oil are given in Table 1. The standard doses were 1.2, 
2 and 3.3 units per day per chick, and the corresponding doses of the 
test material 1.552, 2.587 and 4.312 mg. Thus the dosage steps were 
3:5 and the hypothesized potency of the oil 773 units per gram. The 
response was tarsal-metatarsal distance, after a 5-week experimental 
period. As is usual in this type of test the best metameter was found 
to be log-response, the exact definition being: ‘100 x the 2-figure log 
of the T.M.T. distance, in mm/100, measured on unmagnified X-ray 
photographs of right legs”. It is this metameter that is listed in Table 1. 
For the purpose of illustration, assays and fiducial limits have been 
calculated in full (expression (1) being used for the limits) and the 
latter also obtained on the appropriate chart for (i) all 168 responses: 


TABLE 2 


ESTIMATION OF ERROR FROM THE ASSAY DATA IN TABLE 1 
Percentage Fiduciai 
Selection Mean Limits at P = 0.05. 
of responses n Estimate r A 

(units/g) Computed By chart 
Allitems.... 168 1000 0.150 89.9-112.2 89-113 
| 96 1081 0.158 86.0-118.9 86-119. 
Mext@se .... 48 871 0.139 82.9-122.9 83-123 
. . 24 981 0.149 74.6-143.3 74-145 
Random pairs. 12 1252 0.154 60 .6-284.9 61-280 
Random pairs. 12 954 0.168 54.6-251.7 55-250 
Random pairs. 12 856 0.149 57.7-192.1 57-193 
Random pairs. 12 641 0.127 57.7-153.5 57-154 


(ii) the first 96; (iii) the next 48; (iv) the final 24; and (v) some randomly 
selected pairs of responses from each of the six groups (n = 12). There 
is, of course, a considerable number (in fact, nearly 3 X 10'°) of possible 
groupings of n = 12 from the totality of results; it was thought inter- 
esting to select, by means of tables of random numbers, four different 
ones (i.e. no one response in more than one selection) as examples of 
small assays. The results are juxtaposed in Table 2. It should be 
mentioned that the graphical extraction of the errors was done visually 


. 
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on quite small-scale charts (by an independent observer) without any 
special effort to get “fine” readings. 

To illustrate the degree of approximation attainable by the use of 
the charts for asymmetric assays, the three “sub-assays’” referred to 
above were artificially damaged by the deletion of 7 responses from 
each. The deletions, picked out from tables of random numbers, are 
starred in Table 1. For the smallest sub-assay, as might be expected, 
the asymmetrizing damage is greatest. Consider the medium sub-assay, 
originally with 8 responses per group, and now with 6 + 7+6+8-+4 
7+ 7. To evaluate we begin by calculating the new group averages, 


and the overall standard deviation (now with 35 degrees of freedom). 
We find: 


S, = 142.00; U, = 150.29; S, = 100.00; U, 


107.63; S, = 87.86; 


Us = 94.86;s = 17.75. 


These data are then used to apply the validity tests described in the 
preceding section: it emerges that, of the three, only what we may call 
the pseudo validity test, that for linearity, gives an abnormal result; 
the probability of random deviations’ from a linear response curve 
being so large is apparently less than 0.02. (In the whole assay of 
168 responses, the corresponding probability is between 0.1 and 0.2.) 
We can therefore accept the assay as meaningful. 

All steps to this point are essential, no matter how the error evalua- 
tion is to be tackled. But from here eunreed there are alternative routes. 
To find the exact mean estimate and its limits we evaluate B, V(B), 
D and V(D), with due attention to the weighting necessitated by the 
group inequalities. Then the estimate is calculated in the ordinary 

_ way, and the fiducial limits from 


A(C — 1) + VCA(C — 1) + CPV(D)/B’, 


where A and C have the meanings assigned in the early part of this 


paper. The alternative is to evaluate B and D from the unweighted 
averages: 


= (142.00 + 150.29 — 87.86 — 94.86)/4 = 27.39 
D = (150.29 — 142.00 + 107.63 — 100.00 + 94.86 


— 87.86)/3 = 7.64; 


to calculate the mean potency estimate from B and D: 
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773 X antilog 1D/B = 773 X antilog 0.2219 X 7.64/27.39 


773 X 1.153 = 891 units/g; 


and to make the one further calculation, 


X = sI/B = 17.75 X 0.2219/27.39 = 0.144. 


We then enter the appropriate chart with this \ and n = 41 and read 
off the error range as already described. (‘‘Found/expected potency of 
test material” is here 1.153.) The fiducial limits at P = 0.05 are so 
found to be 81 — 126%. 

A comparison of these alternative treatments for the three asym- 
metrized sub-assays is given in Table 3: 


TABLE 3 
A COMPARISON OF EXACT AND APPROXIMATED EVALUATIONS OF ASYMMETRIC 
ASSAYS. (SUB-ASSAYS FROM TABLE 1 WITH STARRED RESPONSES DELETED) 


Exact Approximated 
n Mean (u/g) P = 0.05 limits Mean (u/g) P = 0.05 limits 
(by chart) 
96-7 1089 85.6-119.7% 1082 85-120% 
48-7 891 81.0-125.9% 891 81-126% 
24-7 913 66 .8-162.8% 921 67-164% 


It will be observed that there is nowhere any appreciable loss sus- 
tained by the replacement of the exact and laborious evaluations by the 
approximate method of arriving at the mean estimate and the use of 
the chart to obtain the fiducial limits. This is not to suggest, of course, 
that gross asymmetry need never be taken into account; on the contrary, 
it clearly must be, the important inference is that the frequently en- 
countered slight asymmetry is biometrically negligible. 


APPENDIX 
I: EFFICIENCY OF ASSAY DESIGNS 


Just as the invariance of an estimate of any quantity can be used 
as a measure of relative value, so the reciprocal of the square of the 
fiducial range of an assay is a measure of its efficiency and has been 
suggested as a routine appendage to assay results. Finney [12] has dis- 
cussed the influence of design on efficiency, and the following treatment 
of 4- and 6-point designs is an adjunct, in slightly different form, to what 
he has demonstrated. 
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Let us confine attention to assays with A = 0, i.e. with accurate 
predictions of potency. The log fiducial range is 


+ IV(C — 1/3; 


so that we may take 
efficiency « 3/I°(C — 1)(k’ — 1). 
This becomes, if C is expanded, and if, in the expansion, s//B is replaced 
by A, 
— — 1) 


It follows that the relative efficiency of the 4-point and the 6-point 
designs for a given n (determinate) and a given \ (a property of the 
experimental animals) is, 


n/t, — 
n/ts — 3d°/2I3’ 


the subscripts designating the appropriate allocations. Note that, for 
the given n, there will be two ?’s because of the difference of 2 degrees 
of freedom associable with the residual standard deviation in the two 
designs [(n — 4) — (n — 6)], a difference whose importance dwindles 
as n grows. The value of this efficiency expression normally depends, 
it is plain to see, not merely on the ratio of the two J’s, but on their 
absolute values; in other words efficiency statements about compared 
designs with a 1:2 dosage step for 6 points and a 1:4 step for 4 points 
will not, strictly speaking, be applicable when the steps are, for instance, 
3:5 and 3:83 respectively; and in general the greater the basic dosage 
step (J;) the greater the efficiency ratio will be. However, for the values 
of n and \ normally dealt with, this influence of J is trivial. And in 
one special case it is wholly absent: when I,/I; = +/8/3, the efficiency 
ratio will be asymptotic to unity (the ¢ difference being the limiting 
factor) whatever the values of J; moreover, in this circumstance alone 
the value of A is without influence on the efficiency. A further point 
to bear in mind is that the efficiency here defined is a function of the 
working probability level, which we shall take as 0.05. 

The tabulation following gives some percentage efficieacies of two 
kinds of 4-point assay relative to a 6-point assay, with J = log 2. The 
first kind, “halved dosage range’’, is a 4-point design with the same /, i.e. 
we can imagine the subjects from one of the end pairs of six dosage 
groups removed and everily re-distributed among the remaining four 


Efficiency ratio, 4-pt: 6-pt = 


‘ : | 
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groups. The second kind, “equal dosage range”, is a 4-point design 
with J = log 4, i.e. as if the subjects from the middle pair of six dosage 
groups were removed and evenly re-distributed among the remaining’ 
four. 


PERCENTAGE EFFICIENCY OF 4-POINT ASSAYS COMPARED WITH 6-POINT ASSAYS 
WHEN THE STANDARD AND TEST DOSES YIELD EQUAL RESPONSES 


Values of n: 

r 10 15 20 25 30 50 100 

0.1 109 96 96 96 97 98 99 

0.12 98 91 93 94 95 96 98 

0.14 83 86 89 91 93 96 98 

halved 0.16 62 78 85 88 90 94 97 
dosage range. | 0.18 32 70 79 84 87 92 96 
0.2 — 59 73 80 84 91 96 

0.3 — — 17 43 56 77 90 

0.5 — — — — — 16 67 

0.1 138 108 104 102 102 101 100 

0.12 143 109 104 103 102 101 100 

0.14 149 111 105 103 102 101 100 

equal 0.16 159 113 106 104 103 101 101 
dosage range.| 0.18 172 114 107 105 104 102 101 
0.2 193 117 109 106 104 102 101 

0.3 —_— 146 121 114 110 105 102 

0.5 — —_— 808 193 149 117 107 


The lower half of the table, referring to ‘equal dosage range”’, is 
the more interesting: we may presume that a switch from 6 to 4 points 
would involve elimination of the middle points and retention of the total 
dosage range. The object of such a switch would be primarily to in- 
crease the amount of information per animal; and it is perhaps remark- 
able to observe that in fact the increase is negligible—except for very 
small and/or very insensitive assays. The increase is greater the poorer 
the “guess” at the potency, i.e. the further A is from zero, but even so 
it is rarely substantial, and as Finney, who has evaluated some effi- 
ciencies for A ~ 0, remarks, “the increase in the fiducial range . . . 
may be a small price to pay for data leading to adequate tests of assay 
validity” [12|—referring to a change from 4 points to 6 points. There 
may however be dosimetric or isogenic reasons in favor of the simpler 
design. 


: 
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II: SYMBOLOGY 


An explanatory reference list follows of symbols employed in this 
paper that are not universally accepted, although most are in some use, 
and several (C, M, R, S, b, g, \) in wide use, in the statistics of bioassay. 
A =D/B 
B = mean response increment between adjacent S’s and between 
adjacent U’s 

C = 1/(1 — g) = Fieller’s statistic. [See ref. 10] 

D = mean response difference between comparable S’s and U’s (as 
U minus S) 

I = log of the ratio of doses given to adjacent S’s or adjacent U’s 

L = aa fiducial limit of A measured from a mean of zero. Thus 
antilog I (A + L) is a fiducial limit (pz or py) of R, while 
antilog (2 + JL) is a percentage fiducial limit of R and of the 
absolute potency estimate of the test material. [See ref. 13]. 

M = IA = D/b 

R = antilog M = best estimate of the unknown potency (=p = 
antilog Ja), in terms of the hypothesized potency of the test 
material = potency ratio of comparable doses of test material 
and standard preparation. 

S = mean response of one dosage group of the standard prepara- 

tion. Distinguished as S, --- S, 
U = mean response of one dosage group of the test material. - Dis- 
tinguished as U, --- U, 

b = B/I = best estimate of the mean linear slope of the two re- 

sponse/log-dose curves 

g = tV(B)/B’ = 1 — 1/C = Finney’s statistic. [See ref. 12] 


k = number of dosage groups in S or U (to form, therefore, a 
2k-point assay design) 

n = total number of responses in the assay 

a = parametric value of A 

= s/b = a coefficient of response sensitivity. 
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Inu HM#emoriam 


Charles P. Winsor 


It is our unhappy duty to report the sudden death of Dr. C. P. 
Winsor on April 4. Charles Winsor was born in 1895 and was educated 
at Harvard University, originally as an engineer. After seven years 
with the New England Telephone and Telegraph Company, his in- 
terests shifted to biological research and biometry, first at Johns Hopkins 
University and later at Harvard where he received his Ph.D. in general 
physiology in 1935. In 1938 he was appointed to the Department of 
Mathematics, Iowa State College and during the war worked in Prince- 
ton as a research associate under a national defense research contract. 
He joined the Department of Biostatistics in the School of Hygiene at 
Johns Hopkins in 1946, a position which he held at the time of his 
death. He was a charter member of the Biometric Society and first 
vice-president of its Eastern North American Region. As editor of 
“Human Biology” and as a biometrical consultant he exerted a wide 
influence upon both biometry and biology. He will be sorely missed 
by his many friends and colleagues. 
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QUERIES 


GerorcE W. SnepEcor, Editor 


QUERY: Nathan Mantel writes as follows: I believe an error 

88 crept into your reply to Query 84 in the September 1950 issue 

of “Biometrics”. In the example you consider at the bottom 

of page 299 you show the following probabilities, pointing out that for 
this sample the exact test is closer to the uncorrected chi-square test: 


Chi-square, P = 0.038 
Corrected chi-square, P = 0.096 
Fisher’s exact test, P = 0.046 


For the two chi-square tests, the alternative hypothesis is merely that 
the mortalities are different; but for the form in which you made the 
exact test, the alternative hypothesis is that the mortality in the treated 
group is greater than that in the untreated. For making the desired 
one-sided comparison, Fisher took advantage of the fact that chi, with 
one degree of freedom, is normally distributed. You can make the 
comparison on the basis’ of a two-sided test by including in the exact 
test all of the more extreme cases; namely, 1, 2, 6 and 7 deaths. . It 
happens in your example that the resulting probability, 0.092, is double 
that for 1 and 2 deaths, making the result almost the same as that got 
from corrected chi-square. 


ANSWER: Thank you for calling attention to this blunder. It is one 
"which I have made before, so I am glad to be set right. 


QUERY: In answering Query 85 (Biometrics 6:435, December 

89 1950), Bliss and Winsor indicated that the range could be used 
instead of the standard deviation. Is this always true? If so, 

we can save a lot of time now used in calculating standard deviations. 


It is always true that the range in a sample can be used 
ANSWER: for estimating o in the normal population from which the 

sample is drawn, but such use would ordinarily be un- 
economical. This subject is discussed by Dixon and Massey in chapters 
15 and 16 of their book, “Introduction to Statistical Analysis.’”” The 
copy I have is an offset edition, but McGraw-Hill has just announced 
publication by their firm. 
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In Query 85, thousands of samples had been taken, each consisting 
of 3 readings on some automobile operator’s blood. For purposes of 
statistical examination, the samples cost nothing. As was pointed out, 
in samples of 3 the efficiency of the range, as compared to the standard 
deviation, is 99%. This means that, for estimating a, 100 ranges give 
the same information as 99 standard deviations. Bliss and Winsor 
observed that for the purposes described it would be a waste of time to 
calculate standard deviations. Similar situations may arise in quality 
control and in establishments where multiple determinations are made 
as routine procedure. 

But in most researches, data are expensive as compared to the cost 
of calculating the standard deviation. Suppose for example that a 
chemical determination costs $2 in time, materials and equipment. 
From tables of the efficiency of the range it is learned that in samples 
of 7 this efficiency is 91%. The meaning is that for estimating o the 
standard deviation in samples of 7 determinations is about as informa- 
tive as is the range in samples of 7/0.91 = 8 (approximately). So, for 
the same amounts of information the calculation of s at a cost of 10 or 
15 cents is equivalent to one additional chemical determination at $2. 

Habitually, I use the range for making rough estimates of ¢. On 
the basis of these estimates I am often able to advise my clients that 
significance or non-significance is obvious. But when decisions are close 
and observations are expensive, the computation of s is ordinarily the 
economical procedure. 
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Members of the Society will be interested in the “Directory of 
Scientific Organizations” recently published in Paris by UNESCO. 
This gives a useful resume of each organization arranged alphabetically 
and starting with The Biometric Society. 

Dr. Lloyd C. Miller, Director of Revision of the U. 8. Pharmacopeia, 
has accepted the Chairmanship of the Committee on Sustaining Mem- 
bers. Those on his committee who will be responsible for inviting 
institutions to Sustaining Membership in ENAR are the following: 
E. B. Carter, A. C. De Graff, M. W. Green, W. E. Hambourger and 
B. L. Oser. 

ENAR: By unanimous vote of the Council of the American Sta- 
tistical Association the Eastern North American Region became the 
first organization to affiliate with the ASA under its new constitution, 
which went into effect in 1948. 

The following committees have been appointed in ENAR: Nomi- 
nating—P. T. Bruyere (Chairman), R. L. Anderson, L. N. Hazel and 
L. C. Miller; Membership—Boyd Harshbarger, Chairman, and 16 mem- 
bers; Regional Advisory Board—J. Tukey, Chairman, and 30 members; 
and Joint Meeting with the Institute of Mathematical Statistics—A. S. 
Householder, Chairman, and 4 members. 

ENAR met jointly with the Institute of Mathematical Statistics at 
Oak Ridge, Tennessee on March 15-17, 1951. In the registration of 
104, 35 were members of the Society. At the morning session on 
Thursday contributed papers were presented by A. E. Brandt, W. J. 
Youden, C. Eisenhart, R. A. Bradley and M. E. Terry, and M. L. 
Clark and F. X. Lynch. The afternoon session considered the use of 
stochastic processes in public health statistics, with papers by E. Fix 
and J. Neyman, and by B. G. Greenberg. The Friday morning session 
concerned bioassay with quantal responses, the speakers being J. 
Berkson, J. W. Tukey and I. Bross. In the afternoon the topic was 
Experimental Design and the speakers were R. C. Bose, O. Kempthorne, 
A. M. Dutton and L. D. Calvin. At the banquet on Friday evening 
E. R. McCrady served as toastmaster, and A. M. Weinberg, Director 
of the Oak Ridge National Laboratory, was the speaker. The closing 
session on Saturday morning presented three papers on multivariate 
analysis by 8. N. Roy, A. W. Kimball and D. V. Tiedeman. Abstracts 
of some of these papers are being published in BIOMETRICS. 

WNAR: The Annual Meeting of the Western North American 
Region of The Biometric Society will be held jointly with the Institute 
of Mathematical Statistics and the Econometric Society at the Rand 
Corporation, Santa Monica, California, June 14 to 16, 1951. 
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Région Frangaise (Société Frangaise de Biométrie): En remplacement 
de Monsieur LaMotte, Monsieur Georges Teissier a été élu membre 
du Conseil Régional, 4 la seance du 7 Mars, 1951. 

A la réunion de la Société le 9 Mai au Laboratoire de Zoologie de 
l’Ecole Normale Supérieure, Paris, |’Ordre du Jour a été “Les possi- 
bilités pratiques de l’Analyse séquentielle’, par Mme. Cassignol et 
M. P. Schutzenberger. 
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